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Abstract 

We extend the notion of generalised Cesaro summation/convergence 
developed in [1] to the more natural setting of what we call "remainder" 
Cesaro summation/convergence and, after illustrating the utility of this 
approach in deriving certain classical results, use it to develop a notion 
of generalised root identities. These extend elementary root identities for 
polynomials both to more general functions and to a family of identities 
parametrised by a complex parameter /x. In so doing they equate one 
expression (the derivative side) which is defined via Fourier theory, with 
another (the root side) which is defined via remainder Cesaro summation. 
For n € Z< these identities are naturally adapted to investigating the 
asymptotic behaviour of the given function and the geometric distribu- 
tion of its roots. For the gamma function we show that T satisfies the 
generalised root identities and use them to constructively deduce Stir- 
ling's theorem. For the Riemann zeta function, the implications of the 
generalised root identities for n — 0,-1 and —2 are explored in detail; 
in the case of /i = — 2 a symmetry of the non-trivial roots is broken and 
allows us to conclude, after detailed computation, that the Riemann hy- 
pothesis must be false. In light of this, some final direct discussion is given 
of areas where the arguments used throughout the paper are deficient in 
rigour and require more detailed justification. The conclusion of section 
1 gives guidance on the most direct route through the paper to the claim 
regarding the Riemann hypothesis. 

1 Introduction 

This paper expands the concept of generalised Cesaro convergence introduced 
in [1] and uses this to introduce a notion of generalised root identities. In turn 
these identities are used to investigate both the asymptotic behaviour and the 
geometric distribution of roots and poles of the gamma function, T, and the 
Riemann zeta function, £. 

In section 2.1 we introduce the core notion of "remainder" Cesaro summa- 
tion/convergence, where the partial-sum function analysed via a Cesaro scheme 
is now on a contour relative to chosen point z , rather than always being on 



[0, oo). As an initial example we show that the Hurewicz zeta function, 

is defined naturally in this context, and that well-known relationships among 

Bernoulli polynomials follow naturally within this remainder Cesaro approach. 

In section 2.2 this is extended in the usual way to a notion of remainder 
Cesaro products and, when applied to the elementary case of f(z) = z, is used 
to give a Cesaro definition of the Gamma function, T, at least after ensuring 
that the geometric location of terms in the remainder sum/product is respected 
in the Cesaro analysis. 

In section 2.3 we demonstrate the utility of this remainder Cesaro approach, 
and in particular its geometric intuitiveness when coupled with the dilation- 
invariance of Cesaro convergence, in analysing three example results. The first 
is an elementary derivation of the duplication formulae for T; the second is the 
derivation of the functional equations for T and £, by exploiting symmetries in 
the Cesaro context to introduce bi-directional sums which are apriori periodic 
and thus amenable to Fourier series analysis; and the third is a short derivation 
of an interesting integral identity for (h ■ 

In section 3 we then turn to the question of root identities, and in partic- 
ular the possibility of generalising the well-known identities relating roots and 
coefficients of polynomials to some more general class of functions. 

In section 3.1 a family of such identities, for arbitrary z € C and \i € 
Z>i, are initially derived on a formal, heuristic basis. In general, these equate 
a quantity defined in terms of the derivatives of a function (the "derivative 
side") with a quantity defined in terms of the roots of the function (the "root 
side"). The latter quantity in fact includes equally contributions from roots 
(with multiplicity M > 1), poles (with M < — 1) and even branch points (with 
e.g. M = i) of a function; all of which cases are thus subsumed thereafter 
under the label of "generalised roots". A natural equivalence is then established 
between satisfying the root identities at a single zq V/i€Z>i, and satisfying 
them at arbitrary zq £ C (at least an open neighbourhood) for /i = 1 alone. 

Unfortunately it is readily seen that most non-polynomial functions do not 
satisfy these generalised root identities. However, at least for those where the 
discrepancy is sufficiently well-behaved (as a function of ^o), we show that it 
is possible to find an "equivalent" function with the same root-set which does 
satisfy these identities. The construction of this equivalent function by iterative 
removal of obstructions to the successive (/i = 1,2,3, . . .) root identities at a 
given fixed z a is outlined. 

In section 3.2, however, we verify that in practice many natural and impor- 
tant functions do already satisfy these identities with no, or minimal, adjustment 
by considering three examples which will become primary focuses of attention in 
this article. In the first we show that f(z) — cos(^) satisfies the identities for 
all z € C and all /i e Z>i and we note that this leads naturally to evaluation 
of some values of (, for example Euler's famous evaluation of ((2) = In the 
second we verify that T(z + 1) also satisfies the root identities for all z G C 
and all /i e Z>i, albeit in the case of \x = 1 that a renormalisation is required 
on the root side to handle uniformly the ln-divergences which arise for any zq 
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and which do not have Cesaro limits. In the third we find that £ alone does 
not satisfy the root-identity for fi = 1, but see by experimentation using the 
first 100, 000 non-trivial zeros that the obstruction in this case appears to be a 
constant, — |ln7r, independent of s . This suggests that instead ir~i((s) does 
satisfy the fj, = 1 root identity for all s and we find that in fact this is true 
and reflects a well-known explicit formula for £ which represents the Hadamard 
product formula for £(s). Thus, in this case, the generalised root identity for 
H = 1 is equivalent to the Hadamard product formula; for fi e Z>2 the extra 
factor 7r~ 2 no longer contributes on either the derivative or root sides and so ( 
itself then satisfies these identities for all s directly. 

In section 3.3 we then turn to the possibility of further generalising our root 
identities by allowing \i to be not just a positive integer, but /i£C arbitrary, 
with the cases of fi € Z<o being of particular interest. We begin by discussing 
extensively how to even make sense of the identities in the case of arbitrary 
complex [i. On the root side a generalised convergence scheme - always Cesaro 
in this paper - is required to make sense of what are now usually divergent sums, 
and we find that it is crucial that the summands, , M i ,„ , in this Cesaro context 
be added in geometrically at the shifted points z — r i- On the derivative side 
the notion of (g^)** now becomes problematic; we note various basic properties 
it should have and then define it consistently with these via Fourier theory, in 
a manner familiar from pseudodifferential operators. 

As a first example of the value of this generalisation to arbitrary complex fj, 
we give a short derivation showing that, for our example of f(z) = cos(^), / 
satisfying the generalised root identities for all z and all fi e C is equivalent 
to the functional equation for (. With this as indication of the utility of this 
extension, we then discuss these generalised identities in the particular cases 
of \x G Z< where the T(fi) factor becomes singular. For the remainder of the 
paper, in fact, we shall ignore the case of arbitrary fi£C and focus exclusively on 
fi G Z< . In this case we show how to interpret our derivative-side definition via 

a distributional interpretation of the expressions r(^+i) ano - n °te that in fact it 
is this singular behaviour that generally makes calculation of the derivative sides 
more tractable for [i G Z<o than for arbitrary /i G C since it allows omission of 
hard-to-compute finite integral contributions, leaving only these distributional 
contributions. We then discuss at some length why, for a function / satisfying 
the generalised root identities, the identities for /i G Z< should naturally give us 
much more information than the corresponding cases of [i G Z>i both about the 
asymptotic behaviour of / and about the asymptotic and geometric distribution 
of its roots. We conclude section 3.3 with a quick overview of the initial formal 
Fourier theory results that we will use in the remainder of the paper, which is 
concerned with analysing the implications of the generalised root identities for 
our remaining two key examples, T(z + 1) and C( s )- 

In section 3.4 we consider the case of T(z + 1). We first verify that T(z + 1) 
does in fact continue to satisfy the root identities for all [i € Z<o by explicit 
calculation of both the root sides (using Cesaro) and derivative sides (using 
Fourier and distributional calculations) in these cases. Then, as an example 
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of the asymptotic information contained within the identities for ji € Z<o we 
show how, for this example of T, we can use them to deduce Stirling's theorem 
describing the asymptotics of T(z + 1). We do this in an engineering fashion 
by successively using the identities for /j, = 0, — 1, —2 ... to construct the next 
order in the Stirling asymptotics. We show, using the generating function which 
defines the Bernoulli numbers, why the Stirling formula in turn leaves the root 
identities for \i e Z>i in tact, without need of further correction terms to remove 
obstructions in these cases. 

Finally, we conclude section 3.4 with some discussion of why our deriva- 
tion of Stirling's theorem in fact implies a "two-sided" result giving Stirling-like 
asymptotics also as z — » — oo, in addition to the actual Stirling asymptotics as 
z — »■ +oo. We use the functional equation for V to verify that this is indeed 
the case, but with specific additional "one-sided" asymptotic contributions as 
z — > — oo, albeit ones which did not derail our derivation above. 

In section 4 we then turn to consideration of the generalised root identities 
for n £ Z<o for (. In section 4.1 we consider the derivative sides of these 
identities and demonstrate that in fact they are all identically zero for all so 
when n G Z< . This is done by using the Euler product formula for £ together 
with a basic property of to in fact derive an explicit formula for the 

derivative side as a sum over primes p. 

In section 4.2, which then takes up the bulk of the remainder of the paper, we 
consider the cases of /i = 0, /i = — 1 and /i = — 2 in turn on the root sides with 
a view to seeing whether these root sides can likewise be seen to be identically 
zero for all Sq in spite of not knowing the exact location of the non-trivial roots; 
and indeed whether any asymptotic or geometric information can be gleaned on 
the location of these non-trivial zeros from the requirement that the root-sides 
be zero in these cases. 

In section 4.2.1 we consider the case of \i = 0, in section 4.2.2 the case of 
[i = —1, and in section 4.2.3 the case of \i = —2. In all cases the evaluation of 
the root-sides proceeds by considering first the trivial roots, T, then the simple 
pole and finally the non-trivial roots, NT. 

As expected, the majority of the effort is in the computations for the contri- 
butions from NT. Here we rely on the explicit Riemann-von Mangoldt formula 
for N(T), the function which counts non-trivial roots with imaginary parts 
between and T; this expresses N(T) as the sum of three pieces N(T) = 
N(T) + S(T) + ±6(T) where N(T) is the leading asymptotic part of N(T), 
S(T) is the famous argument of the Riemann zeta function, and S(T) is a com- 
pletely explicit function giving the negative order asymptotics of N(T). The 
calculations of the root-side contribution from NT always thus further subdi- 
vide into calculations of the Cesaro limits for the terms arising from each of 
these three pieces. In carrying out these Cesaro computations we confirm that 
it is absolutely critical when /x e Z< both (a) to respect the location of the 
roots in determining the geometric placement of summands and (b) to indepen- 
dently carry out computations for the NT-roots above and below the real axis 
using separate parameters T and T (and associated geometric complex variables 
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z and z), but to recombine these into a 2-d partial sum function (of z and z) 
before taking limits. Failure to do so leads to incorrect evaluation and paradoxes 
involving apparent inconsistency between the results for fi = and fi = — 1. 

For fi = in section 4.2.1 we verify that we do indeed have the root-side 
identically zero for all s , at least after assigning Cesaro limit to a residual 2-d 
log-term, In (|) (with this latter seemingly ad hoc calculation either following 
from a suitable 2-d Cesaro analysis, or alternatively to be treated as a conse- 
quence of the root identity for \i = 0, allowing it at least to be used confidently 
for the subsequent cases of /i = — 1 and /i = — 2 where it also arises). 1 

For [i = —1 in section 4.2.2 we likewise verify that we do indeed have the 
root-side identically zero for all sq, at least subject to an estimate on S(T) which 
is known to hold conditional on the Riemann hypothesis (RH). 

For /Lt = — 2 in section 4.2.3, however, we find that a symmetry of the NT- 
roots is broken. For /i = and \i = — 1 the fact that any non-trivial roots off 
the critical line occur in mirror pairs meant that the net contributions from 
any such pairs to our Cesaro calculations was exactly as for a corresponding 
double-root on the critical line, and so no explicit terms expressing potential 
contributions from departures from the critical line arose. But for fi = — 2 there 
is such a term, X e , and after carrying out the analogous, albeit more intricate, 
computations for the root side in this case we find that the value of the root-side 
for fi = — 2 is given by — | + X f . 

In order for the fi = —2 root identity to hold, this would need to be identically 
zero and this allows us to conclude section 4.2.3 by inferring, as a consequence, 
that the Riemann hypothesis must be false. 

Given the nature of this claim, we end by including a final section, section 
4.3, itemising and discussing the admittedly many areas where the arguments 
in the paper up to this point have been informal, heuristic or otherwise in 
need of more rigorous treatment. We nonetheless note that, while our emphasis 
throughout the paper is on developing and applying calculational tools rather 
than on formal rigour, we believe the nature of the resulting computations is 
sufficiently suggestive that we are convinced of their correctness, and certainly 
convinced that they warrant further analysis to try to fill the gaps. 

Finally we note that, for those wishing to focus as rapidly as possible on the 
calculations for £, section 2.3 and the discussion of Stirling's theorem in section 
3.4 may be omitted en route. 

1.1 Notation 

In this paper we use the definition of Bernoulli numbers as given recursively by 



Bq = 1 and 




Vn > 2 



(1) 



1 MVe discuss this issue further in section 4.3 and in [8] we return to address it again and 
justify why assigning this 2d Cesaro limit the value is legitimate 
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which implies B\ — B 2 = g, B4 = Bq = ^ and so on. For n > 1 we 
have B 2n +i = and the even-index Bernoulli numbers are related to the values 
of £ at negative odd integers by i?2n = —2n((—2n + 1). A generating function 
for the B n is given by 



t \ - -B ra 



e 



n=0 

Bernoulli polynomials are then defined by 



J2^ t n (2) 

1 ^ n\ 



B n ( X ) : X (J) BjX n ~ 3 = £ (") 



(3) 



j=o w/ 3=0 

and are monic polynomials satisfying B n (0) = B n = B n (l). Since B n (0) = 
B n (l) we alternatively let B n (x) be the periodic, period- 1 extensions of the 
Bernoulli polynomials from [0, 1] to all of R. For n > 3 these are differentiable 
and satisfy j^B n (x) = nB n -i(x), and this relationship continues to hold for 
n = 2 except at the integer points (of measure 0) where B\ (x) is defined to have 
value 0. The first few B n (x) and B n (x) are given by 



B 1 (x)=x-^ , B 1 (x) = {x}-^ 

1 1 

B 2 {x)=x 2 -x+- , B 2 {x) = {x} 2 -{x} + - 

B 3 (x)=x*- 3 2 x 2 + 1 2 x , B 3 (x) = {x} 3 -l{x} 2 + 1 2 {x} 



(4) 



and so on where {x} = x — Floor(x) is the fractional part of x. 

We denote by b n (k) the polynomials in k e Z >0 given by the sums of (n— l) st 
powers: 

bn(k) := ]T i"" 1 (5) 
i=i 

which turn out to be closely related to the Bernoulli polynomials. 

We take the Fourier transform and inverse Fourier transform to be defined 

by 

m (?) - /_°° fix) ■ e- fa « dx and J"" 1 [ 5 ] (.t) - i- |°° fl (0 • e te « d£ 
Finally we denote the odd Heaviside function by H (x) so that 
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Ho(x) = 




,x > 
,x < 



while Hq(x) denotes the standard Heaviside function 




,x > 
,x < 



2 Remainder Convergence 
2.1 Remainder Cesaro Summation 

In [1] the notion of Cesaro convergence was recast in terms of a Cesaro operator 
P[f](x) := \ Jq f(t) dt (on functions / : R — > C) and extended by consideration 
of the eigenvalues and eigenfunctions of P so that a classically divergent function 
/ has generalised Cesaro limit L if q(P) [/] (x) — > £ as x — >• oo classically for some 
regular polynomial in P, <?(P) (regularity being equivalent to having q(l) = 1). 
Since the eigenfunctions of P, namely x p with eigenvalue (p G C), arise 
naturally in the Euler-McLaurin asymptotic expansion for the partial sums of 
the Riemann zeta series, J2^Li n ~ S i generalised Cesaro summation can be used 
(see [1]) to derive the analytic continuation of ( to Re(s) < 1. 

While various "schemes" were considered in [1], for ( the cleanest extension 
arose from a geometric Cesaro approach in which a partial sum function (of a 
real variable) is obtained by adding in each term n~ s at the point z = n. 

It is more natural, however, to view this analysis in [1] as a special case 
of remainder Cesaro summation in which contributions occur at points located 
relative to an initial point zo € C. Specifically, given a function / : C — > C we 
wish to define its strict remainder sum at z £ C by 

oo 

R + [f](z ):=J2f(zo+n) (6) 

71 = 1 

where we want to understand the sum on the RHS in a generalised Cesaro sense 
along the horizontal contour 7 : t zq + t. Adapting the working definition 
in lemma 4 in [1], we do this by adopting the following working definition of 
Cesaro convergence along a contour: 

Definition 1: Suppose 7 : [0, 00) — > C : t i-> -f(t) is a contour parametrised 
by arc length t, starting at z with j(t) — > 00 as t — > 00, and suppose f is a 
function on 7 which can be written as 

n 
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for some finite collection of constants dj G C, pj G C\{0} and rrij G Z>o and 
some remainder function R(t). If there exists n G Z>o such that P n [R](t) — > L 
as t — »■ oo, then we say that f has generalised Cesaro limit L along 7, where P 
is now the averaging operator (in t) along 7 defined by 



Notes on Definition 1: (i) In the rest of this paper the only contours we shall 
need to consider are either horizontal rays parallel to the real axis (7 : 1 1— > Zo±t) 
or vertical rays parallel to the imaginary axis (7 : 1 1— > zq ± it) , which simplifies 
interpretation. For these cases it is not hard to see uniqueness of limits in 
definition 1. 

(ii) It is critical in this definition, and will be crucial throughout this paper, 
that the functions z Pj (\nz) mj that we "throw away" (i.e. assign generalised Ce- 
saro limit for pj =/= 0) are functions of the geometric variable z — "f(t), not just 
the arc-length parameter t; we continue to think of these as Cesaro eigenfunc- 
tions and generalised eigenfunctions in this setting even though this association 
is now loose given the necessary definition of P in terms of t rather than z. 
Of course, in the case where pj £ Z the need for this distinction disappears, 
at least when 7 is a simple ray as per (i), since in this case a straightforward 
Taylor expansion allows expression of each z Pj (In z) mj as a linear combination 
of functions of the form tP i (\nt) nii with all pi 7^ 0, and vice- versa. When pj 
is an integer, however, this distinction becomes pivotal in order to avoid pure 
powers of In, which correspond to Cesaro "eigenfunctions" with eigenvalue 1 and 
thus have no generalised Cesaro limit (see [1], section 2). 

(iii) In relation to this last remark, recall the key point that Cesaro definition 
of generalised convergence is fundamentally intended as a tool for constructive 
analytic continuation (in zq and some parameter /i G C which drives the values of 
the pj 's in the situation being analysed) . As such, the above distinction could be 
rephrased as saying that when all pj ^ Z the definition 1 can be recast, at least 
for 7 a ray, purely in terms of the parameter t (and then fully justified rigorously 
in terms of genuine eigenfunctions and eigenvalues of a regular polynomial in 
Pt); but when p is such that there is a pj G Z>o we find that we need to 
distinguish between z and t as per the discussion in (ii) in order to obtain the 
correct analytic continuation across this value of p. 

(iv) When 7 is just a horizontal ray then, for Zq real, definition 1 devolves 
simply to the case considered in [1] and so simply represents a natural extension 
of the definition in [1] to achieve the desired constructive analytic continuation 
from 1 to C just mentioned. The case of vertical rays can be viewed similarly. 
While the motivating example of remainder summation in (6) leads to a partial 
sum function with evenly spaced jumps, the function / in definition 1 obviously 
need not have any jumps or even arise from a summation process, or even if it 
does, may arise in a way where the spacing between summands is not even. 

(v) We will return briefly to discuss definition 1 and the above observations 
(i)-(iv) further in the final section of the paper, by which stage the intervening 
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calculations will hopefully have both demonstrated the utility of this definition 
and clarified the meaning and significance of some of these comments (which 
may seem obscure at present). However, since our primary focus in this pa- 
per is on calculation, we defer further discussion here and simply conclude by 
summarising definition 1 into the following "working recipe" for calculation of 
generalised Cesaro limits along contours, as it will be used throughout the paper 
to perform the constructive analytic continuation (in zq and /i) just mentioned: 

(a) first remove linear combinations of Cesaro eigenfunctions and generalised 
eigenfunctions, but critically doing so geometrically as eigenfunctions in z = 
jit), rather than simply t, and then 

(b) apply a suitable power of the Cesaro averaging operator "along the con- 
tour" (i.e. averaging in the contour arc-length t). 

Returning now to the case of remainder summation in (6), here we form the 
contour 7 : t i-> zo+t, form the partial sum function s/(zo, t) :— J2 n <t f( z o + n ) 
and use definition 1 to obtain R + [f](z ) as Climst(z ,t), where z := z a + t 

z— >oo 

and we have used the notation Clim rather than Clim to emphasise that it is 

z— >oo t— >oo 

eigenfunctions in z rather than t that we remove from s/(zo, t) prior to averaging 

as discussed above. 

In similar fashion we define R + ^ [f](z) :— X)^Lo^( z + n )' R-[f]( z ) : ~ 

12™=if( z _ n )> R +,o,-[fK z ) '■= I]r=-oo/( z + n ) and so fortn - In tne last 
case with a bi-directional sum, the definition is via two independent Cesaro 

sums, for i?+,o and i?_ respectively. 

In the case of f(z) = z~ s we see at once that £(s) = R + [z~ s ](0) 2 and the 
analysis referred to in [1] consisted of detailing the generalised Cesaro interpre- 
tation of this remainder summation at for arbitrary s. 3 For z ^ 0, R + {z~ s ](z) 
in turn yields the Hurewicz zeta function (h(z, s) — X)^Li( z + n )~ s vm identical 
Cesaro analysis. 

As a first simple example of the flexibility allowed by letting z € C and mov- 
ing to the remainder Cesaro viewpoint, consider in this context the polynomials 
in k e Z defined in section 1.1 by b n (k) := Y^j=i3 n l ■ These are related by 
formal differentiation of the discrete summation variable, k, namely 

±( bn ( k )) = (n-l)-{b n - 1 (k)-((2-n)} (7) 

from which they may be readily recursively generated and shown to be connected 
to the Bernoulli polynomials B n (x) . This formal differentiation may be rigorised 
in the remainder Cesaro approach by naturally defining b n (z), for arbitrary z, 
as 

z 

b n (z) - ]T.r s := R+[^ s }(0) - R+[2- s ](z) = C(s) - ( H (z,s) (8) 

2 The notation here means, of course, (0) where f(z) = z~ s ; we shall use this nota- 

tional shorthand extensively 

3 Here s is playing the role of the variable fj, in the discussion above; with 20 = held fixed 
the Cesaro convergence is used to analytically continue the expression R+[z~ s ](0) from its 
classical region of convergence when Sft(s) > 1 to all s S C. 
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where s = — (n — 1). Utilising the fact that we may now rigorously differentiate 
with respect to z and that j^(h(z, s) = —s(h{z, s + 1) (after noting that sum- 
mation and differentiation may be commuted in this case in a careful Cesaro 
analysis) we immediately obtain the claimed relationship on setting z — k. 

A second way in which remainder summation immediately expands the scope 
of analysis is that R + (and i?_ etc) are now operators on functions and thus 
amenable to spectral analysis. 

2.2 Remainder Cesaro Products 

The notion of a Cesaro remainder product follows naturally via exponentiation 
of the remainder sum of the logarithm: 

Y[[f}(z):=eMR + [Mf )}(*)) (9) 

R 

For example, for the identity function f(z) = z, we have I1_r[^]( z ) = 
cxp(_R + [ln(z)](z)) and the calculation of i? + [ln(z)](2;) by Cesaro means follows 
readily from the Euler-McLaurin sum formula. Specifically, for any given zq, 
using the notation from [1] where a £ [0, 1) and recalling that z In z and z are 
generalised eigenfunctions of P where z — z + k + a, we have 

k k 

H^a + j) = f ln(z + t) At + * ln(z + k) + C Za + o(l) 



c 
c 



Thus 



(z + k+ ln(z + k) - (z + k) + C Z0 + o(l) 
(z a + k + a) ln(z + k + a) - (z + k + a) + C z , 



Y[[z}(z )=eMC Zo ) (10) 



where 



C Z0 = lim J V ln(z + j) -(z + k+h ln(z + k) + {z + k) \ 

[i=l J 

It turns out this example leads directly to the T function: 



(11) 



Lemma 1: The T function is given directly in terms of remainder Cesaro 
products by 
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or equivalently 



\n(T(z + 1)) = R+[\nz}(0) - R+[lnz]{z) = ^ In(27r) - R+[lnz](z) (13) 

This is readily proven by, for example, comparing Taylor series for the two sides 
of (13), but we omit details here. Note, however, that it is critical in the above 
Cesaro analysis that summands are introduced geometrically at ZQ+j and Cesaro 
eigenfunctions are then viewed w.r.t. this geometric variable z — z + j + a, 
rather than simply the summation index j; this need to respect the geometric 
location of summands in the Cesaro setting will be crucial later in this paper. 

Note that the defining relation for the T function, T(z + 1) = zT(z), is 
self-evident from this definition, as is the fact that T(k + 1) = k\ for k e Z. 
Other basic identities for T (and for other functions) likewise follow directly 
from geometric considerations within the Cesaro viewpoint. 

2.3 Geometric Effects Within the Remainder Cesaro Ap- 
proach 

To conclude section 2 we briefly give three examples to illustrate the way geom- 
etry can be exploited in the remainder Cesaro framework. 

Example 1 - Duplication formulae: The well-known duplication formulae 
for T state that for any iieZ^i 

(2^)^r(z + 1) = n z +ir(— ) • r(— ) • . . . • r(— ) (u) 

n n n 

Using the Cesaro definitions in (12) and (13) this follows almost immediately 
from the dilation invariance of Cesaro summation (see Appendix 5.1). For, 
taking logarithms, (14) is equivalent to 

^1 In(27r) + ln(r(z + 1)) = (z + \) Inn + £ ln(r(^ - <^ + 1)) 

i=i 

which, in light of (13), is equivalent to 

1 n ( — "I 

-R + [\nz](z) = (z+ -)Inn- £i? + [ln^](- - { —^) 

But, by dilation invariance, if we dilate each remainder Cesaro sum on the RHS 
by n (i.e. place each summand ln(^- — + fc) at z — (n — j) + kn) and note 

that ln(- — — r p +kn ) = ln(z — (n — j) + kn) — Inn, then the summands now 
intersperse perfectly (see figure 1 in Appendix 5.3) to yield 

RHS= (z+ -)lnn - R+[lnz - \nn](z) 
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The result then follows on recalling that In n is a constant and therefore R + [In n] (z) 
= (-z-|) Inn (which itself follows immediately on setting s = in (8) and using 
the fact that C(O) = -\ and bi(k) = k). 

It is clear from the structure of the above proof that an analogous duplication 
identity will hold for any function defined by remainder summation R + [f](z) 
(or corresponding product) as long as f(nz) can be related readily to f(z). It 
is trivial for example that 

C*(s, S )=»-«f> (*-(?-Jl, 3 ) (15) 

3 = 1 x ' 

Example 2 - Functional equations: The functional equation for T states 
that 

T(z) ■ r(l - z) = (16) 
sm(7rz) 

or equivalently that 

ln(r(z)) + ln(r(l — z) ) = In 7r — ln(sin(7rz)) (17) 

This can be proven in the remainder Ccsaro framework using (13) by converting 
the pair of remainder sums on the LHS in (17) into a bi-directional sum, R +t0 _, 
and then noting that this is automatically periodic with period 1 and thus 
amenable to Fourier analysis. In more detail, rewriting (13) as 

ln(r(z)) = ^ln(27r)-i? + , [lnz](z), 

consider the diagram in figure 2 in Appendix 5.3, for a given choice of zq. 

Now in the case shown, for Im(zo) > 0, and using — n < 9 < ir as the 
principal branch of In, we have ln(j — zq) = ln(zo — j) — in for all j. Let 
z\ := ((k — zq) + a) be the variable in the partial sum, si, for the Cesaro 
computation of i? +i o[hiz](l — zq) and let z 2 := ((z — k) — a) be the corresponding 
variable in the partial sum, s 2 , for computing R-{\nz](z ), so that we have 



R + o[lnz](l — zo) — Climsi(zi) and i?_[lnz](z ) — Clim si{z-i) 



Then 

k k k 

si{zi) = ^2M.j - z Q ) = ^ln(z -j) -iTT^l 

i=i i=i i=i 

= 52(^2) - ink 

Now we know from (11) that 

1 r 
si(zi) = (-zq + k + -) ln(-z + k)- {-zq + k)+ C- Zo + o(l) - C_ 
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Therefore 

52(^2) = -(^o -k - ^)[ln(z - k) - in] + (z - k) + C- Zo + ink + o(l) 

But 

— 02 In 22 = — (zq — k — a) ln(zo — k — a) 

OL 

= -(zo - k) ■ {ln(z - k) - r-} + aln(z - k) + o(l) 

(Zo - K) 



Thus 



-{zq - k) ln(z - k) + ^ ln(z - fc) + ^ 



c 

52(^2) ~ -2 2 lnz 2 + z 2 + C- Zo + in(z - a) 



c „ . , 1 
2 

1, 1 



C- Za +in(z - -) 



= in(z Q - -) + - ln(27r) - ln(r(l - O )) 
and on adding si and s 2 we have 

i?+,o,-[lnz](z ) - ^(z - + ln(27r) - ln(r(l - z )) - H^M) ( 18 ) 
We see that (17) follows if and only if 

R +fi _[lnz](z) = ln(sin(7rz)) + in(z - *) + ln2 



and on writing 



(_p-ivz \ °° 2-rrinz 

(1 - e 2mz ) = -inz + — - In 2 - ^ 



n 



this in turn follows if and only if 

00 

e 



00 27rinz 



i? + ,o,_[lnzj(*) = -£ — — (19) 



n 

n—l 



Now, by definition, i? +i o,-[lnz](z) is a periodic function with period 1, express- 
ible as a Fourier series 



#+,o,-[ln2](;z) = ^2 a « e 



2-7rm2: 



and so we see that (19) is equivalent to verifying that the Fourier coefficients in 
this expansion, a n , satisfy 

a = {~ 1 n > ^ Z >» 

" lo , uez<„ 
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But 



1 

a n = j R+ fi ,-[\nz]{z) ■ e - 2mz dz 
o 

oo 

= j lnz-e- 2mnz dz 

— oo 

oo 

= j \n\z\ ■ e- 2mnz dz + in J e- 2mnz dz 

— OO -co 

In a Cesaro analysis the second term here gives 

Clim— (l-e- 2 ™* ) = ' neZ *° 

x^co 2mn I , n = 

while the first term, on performing integration by parts and using a Cesaro 
calculation to justify ignoring the boundary terms, gives 



oo 

i r l 

mn J z 



£ -2^nz dz 

2mn 



Now the Fourier transform of the constant function 1 is 27T(5 (£) and multipli- 
cation by z corresponds to i^, so division by z corresponds to — i f d£ and we 
see that this last term integrates to a Heaviside function 



2iri 
2mn 



n e Z <0 
n = 
n e Z >0 



Combining terms, this yields at once the required form for a n and the functional 
equation for T follows as claimed from (19). 

We see that the relationship between values of T at z and 1 — z expressed in 
the functional equation (16) arises naturally from the geometry of the remainder 
Cesaro approach shown in figure 2, together with the Fourier-amenability of the 
resulting periodic bi-directional sum i? +j0 ,-- As such it is again the case that 
an analogous functional equation should hold for any function defined, like T, 
by a remainder Cesaro sum. This is indeed the case, for example, for the 
Hurewicz zeta function (h(z, s). While the Fourier leg of the analysis is difficult 
to perform in this case for arbitrary z, for the particular case of z = it can be 
accomplished relatively directly and the functional equation for the Riemann 
zeta function is then derived, namely 

C(l- S ) = 2 1 - s 7r- s cos(^)r( S )C( S ) (20) 
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For reasons of space, however, we omit details here, particularly in light of 
the fact that we will derive the functional equation for £ in section 3.3 by 
alternative means arising from our introduction of generalised root identities. 
Note, however, that in the derivation using bi-directional sums and Fourier 
theory the connection between values at s and 1 — s arises from the Fourier 
theory rather than from the Cesaro geometry as it did in the example of T. 

Example 3 - An integration relationship for (h and (: In [1] it was 

shown that if x = k + a then we have 

k^co 1 , else 

and that this combined with the Euler-McLaurin sum formula leads to the 
derivation of the values of £ as Cesaro limits of partial sums J2n=i n >> - Since 
k p dk = ( — l) p -=p^, so when ,s G Z< we obtain 

C(s)= / 5>-*dfc 

l l n=l 

and in light of (8) this can be re-expressed as saying that 
o 

J( H (z,s)dz = forseZ< (21) 
-l 

In this last example we show how (21) in fact follows for arbitrary s / 1 
as a consequence again of geometric considerations in the remainder Cesaro 
approach together with dilation-invariance arguments similar to those deployed 
in example 1 for the T-function duplication formulae. 

Lemma 2: We have 

o 

J ( H (z,s)dz = Vs ^ 1 (22) 
-l 

Proof: Suppose initially that Re(s) < 1. Using right Riemann sums, the 
integral on the LHS becomes 

^ n — 1 . n—1 

/ ( H (z,s)dz= lim - VCff (--,*)= I™ - V R+[z- s }(- 3 -) 

_1 3=0 3=0 
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But dilating by n and noting that (-)~ s = n s z~ s and that the summation 
points after dilation all interleave to precisely fill up the positive integers, it 
follows from Cesaro dilation-invariance that this becomes 

lim n'^R+lz-"]^) = lim n 8_1 C(«) = 

n— >oo n— >oo 

This proves (22) for Re(s) < 1. At s = 1 the integral is not Cesaro convergent 
but for all other s the result then follows by analytic continuation. 4 

Having introduced and illustrated applications of remainder Cesaro summa- 
tion, we now turn to the primary focus of this paper, in which it plays a key 
role, namely generalised root identities. 

3 Generalised Root Identities 
3.1 Introduction 

For a polynomial p(z) — a n z n it is trivial that the roots are related to the 
coefficients by J2{rootsr i } r i = — ^p 1 an( ^ similar identities for J2 r T' 771 *= 
Z>i. These identities can be recast in a way potentially applicable to more 
general functions, for example entire functions with Taylor series Y^=i a nZ n , 
as relations involving the reciprocals of the roots as follows: Suppose f(z)= 
Y\(z — ri) Mt , then ln(/(z)) = ^ Mj ln(z — rC) and therefore in general for any 
jU € Z >0 and any z 

In this formula note that nothing requires Mi € Z>o- The multiplicities, Mi, 
may be negative integers (n a pole) or indeed arbitrary complex numbers (r» a 
branch point). We shall continue to use the term "roots", or sometimes "gener- 
alised roots", to cover all these possibilities. 

Note also that while (23) gives a family of identities for \i = 1, 2, 3, ... at 
any given z , this is of course equivalent to knowing the single identity for fj, = 1 
at general zq. For if we have (23) for all \i e Z>i at z then for n = 1 at any 

4 Note that the reasoning here likewise immediately yields a family of special values of 
for example, for any prime p, "}2Ch{— ^,s) = (p s ~ 1)C( S ) where the sum is over 1 < j < p; 
for p, q prime, Yl C-tf ( — s ) = ((Pl) s — p s — q s + 1)C( S ) where the sum is over 1 < j < pq, j 
coprime to pq, and so on. 
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z + h in (23) we have 

LHS = - A(i n(/(z))) \ z=Zo+h = - ^(M/»)) | Z=Z0 where /» = /(* + ft) 



1 -dz ( e ^( ln (/W)))U=,o 

- -E^(£fW)>)i. 



J=0 {roots r;} 



, , ( z o - n) 

{roots ri) 1 — 



= - V* — = RH S 

{roots ri) 

and vice-versa. 

Unfortunately, of course, while both sides of (23) now make sense for more 
general functions, the identity is generally untrue for non-polynomials. For 
example, for f(z) = e z , / has no roots so the RHS of (23) is zero while the LHS 
is identically 1 at /j, = 1 for any zq. 

However, suppose / fails (23) at fi = 1 with error function 

{roots ri of f) 

and g entire. Then if we set h(z) — e G ^ f(z) where G'(z) = g{z) it follows from 
the fact that the root-sets of / and h are identical that in (23) for h(z) we have 

LHS = -±MHz)) = -±(G(z) + Mf(z))) 

= —g{z) — — = - > = RHS 

f(z) . th Az- n ) 

{roots ri of h) 

Thus, even though / fails the generalised root identity (23), there is a unique 
(up to an overall scalar), nowhere-zero entire function e G ^ whose product with 
f(z) yields a function with the same root-set satisfying the identity. 

In terms of equivalence classes, if we set two functions equivalent, / ~ h, if 
and only if there exists an entire, nowhere-zero function, k(z), with fc(0) = 1 
such that h(z) = k(z)f(z), then within any equivalence class there is a unique 
representative which satisfies the generalised root identities (23). 

This can be viewed in a different way by thinking of zq fixed and ft = 
1, 2, 3, . . . successively. Taking z a = for example, if the LHS and RHS of (23) 
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differ by a\ for fi = 1 then we can remove this "obstruction" by multiplying / 
by e aiZ since this leaves the root side undisturbed and contributes the required 
ai to the derivative side. Similarly, if the obstruction of the /j, = n identity at 
z = is a n then multiplying / by exp( a "^ ) again leaves the root side of (23) 
unchanged but contributes precisely the required correction of a n to the n = n 
identity, while contributing nothing to any of the other identities with [i e Z^ n 
and thus leaving these identities undisturbed. In this way, working successively 
through fi = 1, 2, 3, . . . we can correct each obstruction and produce a new 

2 3 

function h(z) = cxp(aiz+^|^ + ^|^ + . . -)-f(z) which has the same generalised 
root-set as / and does satisfy the root identities (23) for all ji e Z>i. 

Note that if the error function, g(z), obstructing the fi = 1 identity for / is 
not entire, then the situation becomes more complex. 



3.2 Some Examples 

In reality the root identities (23) hold without any need for adjustment for many 
well-known functions. We next briefly consider three examples. 



Example 1 - /(z)=cos(^): Here the roots of / are all simple roots of mul- 
tiplicity 1 at the points ±(2fc — 1), /c e Z >0 , and so for fi = 1 at arbitrary z in 
(23), we have 

°° 2? 
RHS = -E sg-(2*°-l) a - ^ 

LHS = |tan(^) 

These are equal by a well-known identity, and so f(z) = cos(^) satisfies the 
generalised root identities (23) for arbitrary zo and all /it £ Z>i. 

Since the roots of / are closely related to the positive integers, the RHS of 
(23) leads easily to values of ((n), n e Z >0 . In particular, for \i = 2 the root 
identity at z = immediately captures Euler's famous formula that X^^Li ^ = 

2 

This is an exampe of an "inverse" use of the generalised root identities, where 
a function of interest, £, arises naturally from the root side of the identities for 
a known function, f(z) = cos(^), and information about £ is then deduced by 
considering the derivative side of the identities for /. We shall return to this 
example, along these lines, in the next section. 



Example 2 - T(z + 1): The roots of T(z + 1) are all simple poles (Mi = -1 Vi) 
at z = — 1, —2, —3, Thus for /i = 1 and arbitrary z in (23), we have 

rhs = y -J— 

This is divergent for all z and is not even Cesaro summable because the partial 
sums involve the generalised Cesaro eigenfunction with eigenvalue 1, lnz. How- 
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ever, since this Cesaro-obstruction is uniform, it can be "renormalised" away 
uniformly as 



^ = E(^:-r)+7 (24) 

n=l 



1 1 



^ \ zo+n 



where 7 s=s 0.577 is Euler's constant and arises because the difference between 
the In z divergence we are removing and the partial sums of ^2 „ approaches 7 
in the limit ( lim (^2^=1 ~ ~ m -^0 = 7). 

N— >oo 

But (24) is in turn a well-known expression for — w^TY U=z anc ^ thus we 
see that, after renormalisation of the In Cesaro divergence, + 1) does satisfy 
the generalised root identity (23) for p, = 1 and arbitrary zo (note that the 
expression in (24) retains the required simple poles at Z <0 ). 

The need to perform this renormalisation adjustment arises only for the case 
of p = 1. For p = 2 the RHS of (23) is — X^^Li (z +n) 2 ' w hi°h i s classically 
convergent for all z a and clearly equals the the LHS, namely — ^2 (ln(r(z + 

l)))|z=z ; on differentiating the expression for — r^+ij* U=zo m (24); likewise for 
/' = 3,4 

Thus overall T(z + 1) satisfies the generalised root identities (23) for arbitrary 
zq and all p £ Z>i, albeit after requiring an additional renormalisation when 
p = 1 to uniformly remove the Cesaro non-amenable ln-divergences in this case. 
It is easy to see that the same continues to hold true for the case of general 
T(az + b). 



Example 3 - ((s): The root set of ((s) consists of the trivial zeros at s = 
—2, —4, . . . (with Mi — 1), the simple pole at s = 1 (with Mj = —1) and the 
famous non-trivial zeros in the critical strip < Re(s) < 1 (for which Mj > 
are unknown in general). We let T denote the set of trivial roots and NT the 
set of non-trivial roots, which we will also denote by pi rather than rj. If the 
Riemann hypothesis is true, then the pi occur in conjugate pairs solely on the 
critical line s — \ ; if not then some occur in quadruples via reflection in the real 
axis and critical line: p il pi, 1 — pj, and 1 — pj. 

Of course, the roots in NT are not known exactly and thus tackling the 
root side of the generalised root identities (23) directly is difficult. However, 
working empirically first, we can take a list of, say, the first 100,000 non-trivial 
zeros and use them to test experimentally whether £ seems to satisfy (23) for 
fi = 1. In doing this, however, a renormalisation analogous to the last example 
will have to be carried out to handle the trivial zeros, T; considering the case 
so = initially, this is equivalent to setting ^ T (q~t) = \ (the \ factor arising 
naturally since T only covers the negative even integers). Thus, taking initially 
only the truncated non-trivial root-set NT given, for example, by Odlyzko in [2] 
and consisting of the first 100,000 non-trivial zeros (actually 200,000 including 
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conjugates) we find that the root side of (23) for fi = 1 at s = becomes 
RHS « -C- + 1 + g T^y^yl « -1-2655342 

ATT 

(on performing the last computation numerically to obtain J2nt p~ ~ 0.0230737). 

By contrast the derivative side in (23) gives = - ln(27r) w -1.8378771. 

At once we see that £ does not directly satisfy the generalised root identities 
for fi = 1 at So = 0, with the obstruction in this instance being, numerically, 
0.5723429 w 0.5723649 = § Ihtt. 

Next consider so = 5. Here, by the symmetry outlined above, X^att ( 1 _ p .) = 

and so RHS = — {X)t 7TZ — t — ti _ii }■ I n thi s case ^ ne calculation of the 
renormalised value of the first sum is best done by noting that the trivial zeros 
occur at the same locations as the poles of T(| + 1) and thus, up to an overall 
factor of -1, this sum can be evaluated using the generalised root identity for 
r(f + 1) at s = I namely as |^ = - +7>- We thus 

obtain RHS = -\-\~\ In 2 while LHS = = 2 f i ln2 l ln ( 7r )> 

so that again we obtain the same value of | ln(7r) for the obstruction. 

This suggests, along the lines discussed before, that while ( does not directly 
satisfy the generalised root identities (23) for /i = 1, the function 7r~i£(s) may 
well do so. 

This conjecture can be expressed in a different way. Recall that the func- 
tional equation for ( given in (20) can be re-expressed as simply 

S{s)=Z(l-8) where £(s) := (1 - s)r(| + 1)*"* C(s) (25) 

Since the generalised root identities (23) hold for (1 — s) and T(| + 1) (as dis- 
cussed in example 2), and since a product of functions satisfying these identities 
will also satisfy these identities, so the conjecture that Tr~i((s) satisfies them 
for /j, = 1 and arbitrary s is equivalent to having £ satisfy (23) for /j, = 1 and 
arbitrary s . 

In fact this conjecture turns out to be true. In light of the discussion of T 
in the previous example it is equivalent to the following corollary which can be 
found in [4, pg35]: 



Theorem 1: Let NT denote the set of non-trivial zeros of Q, and NT + denote 
the subset with imaginary part > 0. Then, with sums understood to include 
multiplicities, we have 



ill, 1 r'(f) 

(2s 1) Y,p£NT + (g_p)( s _(l_p)) 2 m7F 
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Thus, overall, both 7r~ 4 C( s ) an d £,(s) satisfy the generalised root identities 
(23) for /i = 1 and arbitrary so, and hence also for arbitrary /i £ Z>i. Since 
the factor only contributes to the root identities when /i = 1, note that for 
£ Z>i the root identities will in fact be satisfied directly by ( (which can of 
course be readily spot-checked experimentally). 

Note that, as discussed in [4], this corollary is in fact the expression of the 
Hadamard product formula for £ (which is a holomorphic integral function of 
order 1). We thus see that the generalised root identity (23) for /x = 1 is in this 
case equivalent to the Hadamard product formula for £. Since the identities 
for /i £ Z>i are just derivatives of the /i — 1 identity, they will contain no 
additional information, but we turn now to considering cases of more general 
fj,, in particular \i £ Z< , in the hope that they will lead to new tools with 
additional content beyond the Hadamard identity. 



3.3 Further Generalisation of the Root Identities 

In (23) our current root identities constitute a generalisation from the case of 
polynomials, relating derivatives of ln(/(z)) at z to sums of integer powers 
of shifted reciprocals of (generalised) roots of /. It is natural, however, to ask 
whether they can be further generalised by allowing /i to be an arbitrary complex 
number rather than just a positive integer; that is 

(ln(/(,)))U =Z0 =e^ ]T j^Ty, (27) 

{roots ri } 

The LHS here we call the derivative side of the root identity and denote 
by df(zo,fi) (or just d(zQ,n) if the context is clear). The RHS is the root 
side and is denoted by r/(zo,/i) (or just r(zo,/i)), and (27) is interpreted as 
asserting identity of these two functions of two complex variables. Of course, in 
attempting such an extension the question immediately arises of how to interpret 
either side of (27): 

(a)(i) To overcome the potential divergence of the sum on the root side (e.g. 
when Re(fi) becomes negative) the RHS must be interpreted via a generalised 
convergence scheme, which analytically continues the RHS from its region of 
convergence in the /x-plane. In all cases in this paper this will be a generalised 
Cesaro scheme. 

(a) (ii) In implementing this Cesaro approach, the sum on the root-side will 
need to be interpreted geometrically. Specifically, each term ^J^^-y mus t be 
added in at the shifted point z — itself in the complex plane (rather than, 
for example, always being added in at itself). Thus, as zq varies, not only do 
the summands vary, but their locations move too. We shall emphasise this by 
writing E {zo -roots r i} on the RHS in (27) ■ 

(b) (i) On the derivative side in (27), the interpretation of ( also becomes 
problematic. Two facts that should hold in any such definition, however, are 



TU) \dz 
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that 

dz 



{a z )\ z=za =a z °{\nar (28) 



and 



, , ^)U=o= ' ' " (29) 

/ 10 , else 

We shall formalise the definition of the LHS of (27) in a manner consistent with 
these criteria below. First, however, we demonstrate the extra power obtained 
from this extension of the generalised root identities to arbitrary complex /i 
by showing how the full functional equation for (, given in (20), follows very 
simply from applying them to the function f(z) — cos(^) considered before in 
example 1 in section 2.1. For in this case, taking zo = in (27) we have 

RHS = r/(0, m) = e raAI (l + e- <ir ")(l - 2"")C(/u) 
while for \x ^ 1 we have 

LHS = d f (0,n) = 




| 2=0 



e 2 7r p 

' (l-2")C(l-,x) 



r(M) 

so that setting df(0,(i) = r/(0,/x) we obtain 



C(i-^) = 2 1 -^cos(^)r( M )C(^ 

which is precisely (20). Thus the functional equation for ( is equivalent to 
cos(^f) satisfying the generalised root identities (27) (after also extending easily 
to the case of \i — 1). 

With this example as a motivation, and (28) and (29) as guidance, we for- 
malise the meaning of the complex derivative in the LHS of (27) via 
Fourier theory, in a manner familiar from linear partial differential operators, 
as follows: 

(b)(ii) Writing / as the inverse Fourier transform of its Fourier transform, 
namely 



g(z) = ^ff°° ' ' Nb-di- 



22 



we define (g(z))\ z=Zo by 

= 2nj J^^KOe" ^ (30) 

With (a)(i), (a)(ii) and (b)(ii) the generalised root identities (27) now have 
a well-defined meaning for arbitrary complex /i. We have considered their in- 
terpretation for f(z) — cos(^) and will examine them next for T(z + 1) and 
C(s), but first we make several general observations about these identities and 
their potential uses: 

(i) While the cases of /i e Z >0 give information on the sums of integer powers 
of shifted reciprocals of roots, (27) gives a much more sensitive relationship 
between the distribution of these roots and the behaviour of the derivatives of 
the log of the function. 

(ii) In paricular, for ^ = 0, (27) should give information about the Cesaro 
count of the roots (adding M, at each point (z — r,)), while for \x € Z <0 we get 
information regarding the Cesaro sums of first, second and higher powers of these 
roots. In particular, taking Re{ji) sufficiently negative should give information 
about the asymptotic distribution of the roots and the cases \x e Z< will be of 
particular interest. 

(iii) For /i £ Z<o, however, it is immediately clear that further care will need 
to be taken in (27) owing to the poles of T in the factor on the LHS. At first 
glance these appear to make the df(zo,/j,) always zero whenever /i G Z<o which 
seems problematic (e.g. in interpreting the case \i — as yielding a count of 
roots). Thus, in general, we will need to interpret (27) distributionally, and for 
fi € Z<o the distributional result we shall rely on in this paper is the following 
fact (see Hi'^cermander's treatment in [3, pp74]) regarding the function 

x a = { Xa ' X> ° 
+ ' [0 , x < 

and its normalised counterpart \+ '■= r(a+i) x + : 

Lemma 3: For n G Z> we have 

X? = 6<?- 1 \x) (31) 

where So is the usual delta distribution. 

(iv) The fact that (27) naturally gives (Cesaro) asymptotic information 
about the distribution of the roots of / makes it a natural tool with which 
to investigate this asymptotic behaviour. Another way of looking at this is to 
compare with the case of (23) in which fj, was restricted to Z> . From an exper- 
imental point of view, since for any positive integer \i the RHS of (23) or (27) 
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diverges as zq approaches any finite root , so the cases of such positive integer 
fj, give a tool for detecting the location of finite roots by looking at the domi- 
nant behaviour of these identities as we let z approach such a location. They 
give little immediate guidance, however, about how the root locations behave 
asymptotically as they tend to oo. 

By contrast, for fi = the RHS in (27) is clearly insensitive to the location 
of finite roots, but its Cesaro calculation depends critically on the asymptotic 
location of the roots as they approach oo; and the same is qualitatively true for 
/i = -1,-2,-3,.... 

In part this in turn reflects the non- locality of the derivatives ( j-) on the 
LHS whenever (i ^ 2^>o, which can briefly be seen by considering an alternative 
approach to defining such derivatives, namely: let Th be the operator of trans- 
lation by h (i.e. Th[g](z) = g(z + h)). Since ^ = lirn(^)(l —Th), so in general 

formally 

lim {e i7r "/i-"(l - TO") 

ft (;)«-...)} 

It is only for \i G Z>i that this series definition truncates and leaves a purely 
local resulting function; for all other fj, the resulting function, (^) (g(z)), 
depends on g(z) as z — > oo and so is intrinsically non-local. 5 

(v) In using the Fourier definition of the LHS in (27) we shall use the standard 
Fourier relationships repeatedly, namely 

?[xf(x)](o = tr[/]'(o & WW]® = mm) 




F[f(x)e iax }(0 = FIM -a) & T[f{x + o)](0 = T[f](0^ 
Since = 2nS (^) it follows that, respecting oddness/evenness, we have 



and 

and in general 
T 



= -2mH (Z) = -2™ ■ {H+(0 \) 



= -2nH (0 ■ £ = -2n(H+(0 \) ■ £ 



5 For example for fj, = — 1 we get (g^) = ^ im { — ^(1 + T h + 1% + ■ ■ -)}which in the limit 
yields the left Riemann sum definition of — / z °° g(u)du. 
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Also, inverting the relationship for Fourier transform of a derivative and apply- 
ing it to the formula for above, we have 

nnx] = _ 27r Ml = _ 2 ^ g o + (0-|) (33) 

and 6 

XT i 1 o • o • - j) p /„.-. 

Jja; In x - x\ = 2-ki — = 2m 12 & ■ • ■ (34) 

with corresponding results for higher anti-derivatives on the LHS. It is these 
latter identities that will be used most directly in what follows. 

From the point of view of calculation on the derivative side of the root 
identities (27), note that Hq(£) ■ £ a = so that, for a G Z <0 , the first term 
in the brackets in these expressions leads naturally to use of the Hi'^cermander 
result (31). The second half of these expressions, , leads, under inverse 
Fourier transform, either to ((-function type contributions (for a G Z> ) which 
may be ignored for zq ^ 0, or to Heaviside-type contributions in zq (for a G Z<o) 
whose finite values for any z cancel to under the factor on the derivative 
side when \i g Z< . As such, for our calculations when \x G Z< in the root 
identities for T in the next section we may effectively ignore these second-half 
contributions and elide the distinction between Hq(0 ' an d £+• 

We now take the two examples of T(z+1) and C( s ); considered in the previous 
section for fi G Z>o, and reconsider them for arbitrary ji G Z<o- We devote the 
rest of this section to the consideration of T, and then turn in the next section 
to focus exclusively on the case of (. 

3.4 Example 2 - The Case of V{z + 1): 

For T(z + 1), on the root side of (27) we clearly have at once 

rr(z ,M)-e^ ]T - Mi = -e""Cg (z„, aO (35) 

{zq — roots ri\ 

When zq = this reduces to the zeta function 

r r (0, M ) = -e^C(M) (36) 

and since (h(zo, /•*) = C(m) — J2j=i i - ^ so f° r example at = 0, — 1, —2, —3, . . . 
we get 

technically in fact .F[ma:] = — 2tt — 2ir~/8o(£) with the extra term arising from the 

constant of integration; however we ignore the additional <5o(£) term here since it will not 
contribute to the derivative side of the root identities at fi £ Z<q in any of our remaining 
calculations in this paper on account of the factor whose denominator diverges at such 

/i. 
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1 2 1 



r r (z Q ,0) = z + 
r r (z , -1) 

rr(*>,-2) 
r r (z ,-3) 
and in general for n <G Z>o 

r r (z ,-n) = (-l) n {b n+1 (z ) - C(-n)} 



1 

12 



1 3 1 2 1 

3 Z ° + 2 Z ° + 6 Zo 

_I 4_ 1 3_ 1 2 

4 Z ° 2 Z ° 4 Z ° 



1 

120 



(37) 



n + 



-B n+1 (-zo) (38) 



On the derivative side, noting (jj)^ = (37)^ 1 (37) an d recalling the identity 
in (24) we have 



1 1 

"27rGu) 



(ln(r(z + l)))| z=Z0 
3 ( 00 

rcr 1 -E 



n=l 



1 



n + .t n 



For all except /i = 1 we can omit the do (£ ) terms arising from the - and 7 terms 
here and thus, by the Fourier identities canvassed in (v), we obtain 

dri "-i' ) = r^/j i «)""'|E- i *»K) e '" £ } e '"° Sd « 

on ignoring the contributions from the | term in J?o(0 = ^o~(0 ~ § as j us t 
discussed. Thus 



^r(2 ,M) 



M-2 el ( ZQ + l)C 



2! 3! 



(ii) + (Ml! 



2! 



3! 



+ ... 



d£ 



f l + i[(^o + i)]e+[-^o 2 -^o-^]^ 1 



/>OG 



+ [^I Z + T2 Z + 2^ Z + 127)] + • • • 



2G 



Now suppose /j, e Z< . Working distributionally as discussed earlier in (iii), 

since for \i e Z< we have fT^rpj = (0 (by (31)); so for M = we have 

immediately LHS Zo fi = 20 + 5, m agreement with the first result from (37). 
Similarly, for \i = —1,-2 and —3, after disentangling the impact of the 
factor and the factorial arising from differentiation, we get agreement with the 
formulae in (37) (each arising as a scalar multiple of the appropriate coefficient 
of £+ -At in the bracket) ; and in general 

V/i e Z< d r (z Q ,n) = rr{z ,n) 

so that T(z + 1) does satisfy the generalised root identities (27) for all fi e Z< . 

Thus, overall, T(z + 1) satisfies the generalised root identities (27) for arbi- 
trary zq and fj, €Z; and in fact, as noted for /1 € Z >0 in a previous section, it is 
easy to adapt the above arguments to verify that this remains true for T(az + b) 
for any a, b eC. 

To conclude our analysis of T(z + 1) we illustrate the power of these new root 
identities for [m e Z< to give information about asymptotic behaviour. We do 
this by constructively deducing Stirling's theorem from successive consideration 
of the identities for /j, = 0,-1,-2,. . ., namely: 

Stirling's Theorem: As z — > +00 we have 

T(z + l) = V2^z z+ h- z e J ^ (39) 
where J(z) has the asymptotic expansion 



00 



T( \ — V* ^ 2n ^ ^ ^ ^ ^ 

(Z) ~ Z** (2n - 1) • 2n^T ~ u z ~ 360^ ( > 

n=l x ' 

Equivalently, as z — > +00, we have asymptotically 

ln(T(z + l)) = (z+^)lnz-z + J(z) (41) 

To deduce this consider first the generalised root identity for [i = 0. By (37) 
we have seen (from a Cesaro count of roots on the RHS) that this means 

As we move through /j, = — 1, — 2, . . . the corresponding polynomials in (37) are 
successive integrals of this (as one would expect from each additional application 

°f (<T~) ' ^' e ' integration, on the LHS), with the only new information each time 
being the value of the integration constant, namely (— 1) M_1 C(^) (by (38)). Our 
requirement is thus to understand (a) what factors must occur in ln(r(z + 1)) in 
order for the terms zo and \ to arise on the LHS in (42) when fi = 0, and then (b) 
how the successive integration constants just identified can be made to appear 
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on the derivative side of the root identities (27) when /i = —1,-2,... by the 
inclusion of further terms without disturbing the identities already examined. 

Now considering the definition in (30) we see that integer powers of z 
will arise when the action of ^ f^° 00 (iQ' J 'F[g](£) • c?£ on e lz °^ consists of dif- 
ferentiation to some integer order. In light of (31) this means we require 

e^sKO = x+ n (0 = r ( i-n) g q } for some n e z >o- We thus re q uire -HsKO 

itself to be of the form C H ^ for some suitable C and p (ignoring the contri- 
bution from the \ in £To(£) = Ho(€) — \ as usual). In particular it is clear that 
when fj, = 0, in order to get a constant term \ we need to have T\g\ (£) = \ 
(so that we get \x+~ > l - e - \&o{€ii acting on e" ^); while to get a term z 
we need to have .F[<7](£) = (so that we get x+ 2 ; i- e - ~^o(£)j acting on 

e tz °£). But, comparing with (32), it is clear that this means we need to take 
g(x) = (xlnx — x) + \ In a; and so (z In z — z) + | In z must appear in ln(r(z + l)). 

These in turn then generate the terms — — \zq in the /i = —1 identity in 
(37), leaving only the ^ term to account for. For this, in the same way we used 
terms e az " to "heal" obstructions for root identities with [i e Z >0 , so we now 
naturally consider terms of the form e^"" to "heal" corresponding obstructions 
when fi G Z< . Taking e? first - on taking logs this becomes g{z) = - and, by 
(32), F[g]{£) = -2maH {£). Thus - when A* = -1, = (eliding 

#o(£) and Hq(£) as usual) and on writing = p^rn we see fr° m (31) that 
we will get a contribution of —a from this term in the root identity (arising from 
— a#o(£) acting on e" ^). In order to match the required ^ we thus take a = 
and so include an extra factor of eA; in T(z + 1) (or yg - in ln(r(z + 1))) in 
order to satisfy the /i = — 1 root identity. Importantly, note that this extra term 
does not disturb the existing identity at p — since there it leads to a finite 
Cesaro integral which is cancelled to by the j^gy factor in the \i = identity. 

So far we thus have contributions z z +2e~ z eAi in T(z + 1), thereby ensuring 
the /i = and /j, = — 1 identities are satisfied. In the same way, working 
inductively, for fj, = —n the existing terms generated up to ^ = — n + 1 will 
integrate to yield the polynomial (— l)"6„ +1 (z ) in the expression (38) for the 
[i = —n root identity (after adjusting for the change in ^jj/j factor on the 
derivative side) and we need to introduce a further term e^ in order to match 
the new constant term — (— l) n ((— n) which is all that remains unaccounted 

for. But then taking logs yields g{z) = with T\g](t) = 27rae~^" "^S' 1 
and so - ^ ± (ityFlg] (0 = -a ■ f^)^ = -anS (Q after simplifying the 
two r terms by cancellation (or alternatively using the functional identity (16)) 
and retaining only the Hq (£) term from Ho(£) as usual. Acting on e zz °^ we 
thus need to take 

a= ™ in order to match the required constant term. 
When n is even (i.e. fi = —2, —4, . . .) we are at a trivial zero of £ and a = 0, 
so that there are no even order reciprocal terms in the asymptotic expansion 
for J(z); when n = 2k — 1 is odd we know that ((—(2k — 1)) = — so that 
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a = 2 fc.(2fc-i) ano - we nave to include a term e 2k t 26 - 1 ) * 2fc -! in our expansion 
for r(z + 1). And, again, each new term included does not disturb our existing 
identities for fx = —n + 1, — n + 2, . . . , — 1, since there the e^ term leads to 
finite Cesaro integrals on the derivative side which cancel to against the p^y 
factor in the fx — j identity, — n + 1 < j < 0. 

Consideration of all the root identities for fx e Z< thus leads us construc- 
tively to deduce 

T(z + 1) = C z z+ h~ z e J ^ 

exactly as in (39) and (40), and it remains only to show that C = \/2n. But, 
letting z = k — > oo through the positive integers, taking logs and noting that 
J(k) — »■ as k — > oo we get 



lnC= lim Un(kl) 

k— >ca 



(k+ ^)lnfc- k 




(k + In A; - k 



But this last expression was evaluated in [1] (in the course of calculating C'(0)) 
as precisely \ ln(27r) and so C — V2tt as required. 

This completes the constructive derivation of the expression in (39) purely 
from the root identities for fx G Z<o satisfied by T(z + 1). To complete the 
proof of Stirling's theorem it remains only to verify that no further terms of the 
form e az , n G Z>i, are required in order to also satisfy the root identities for 
fx E Z>i, i.e. that the expression \/2tt z z +h e ~ z e J ( z ) already satisfies all these 
root identities for T(z + 1) without requiring the healing of any obstructions. 

For ji = 1 we immediately encounter a problem from the term z z e~~ z ; on the 
derivative side of the root identity, after the usual passage from Hq to Hq , this 

yields an integrand of i • = — ~|~ an d this leads to a Cesaro divergent 
integral (because log-divergent) in its action on e lz °^ since we no longer have 
the coefficient to normalise this into a 6- function. However, this merely 

mirrors the corresponding Cesaro divergence of i^qrj on the root side, for 
which we needed to perform the renormalisation described earlier in verifying 
the fx = 1 root identity; a corresponding renormalisation is of course required 
on the derivative side. 

With this in mind we proceed with the other terms: z 2 gives integrand con- 

B 2fc i 

tribution — (£),and each term e 2k ■< 2fc - 1 ) ^ 2k ~ 1 gives integrand contribution 

2F§tT) ' (-*) 2fe_1 (£i(^) = -< 2fc ^f^+ _1 ( b y (32))- Thus, for fx = 1, on 
the derivative side of the root identity we get 



1 / 



1 1 .A B 
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Recalling the generating function for the Bernoulli numbers in (2) we see that 
this becomes 



where the last step follows from a direct Cesaro computation of the integrals. 
Since this formally equals the root side of the /i = 1 root identity we see that, 
(up to renormalisation on both sides), v / 2~7rz z+ 2 e~ z e J ( z ' does indeed satisfy the 
root identity for T(z + 1) without requiring any further correction terms of the 
form e az . 

For [i = 2 and higher, no renormalisation issues arise. At \i = 2, on the 
derivative side, z z e~ z contributes — , z 2 contributes W, e~ z makes no con- 

B 2k 1 , 

tribution and e 2k < 2k ~ 1 '> ^ k ~ 1 contributes — B 2 k am after 2fc-fold Cesaro inte- 

gration by parts. We thus end up with ^ + ^ — J2m=i J^" 1 which is a 

well-known expression for — 4^ (ln(T(z + 1))) \ Z — ZQ . Since we have already ver- 
ified that this equals the root side of the generalised root identity for T{z + 1) 
(i.e. that T(z + 1) does satisfy the /j, = 2 root identity) so it follows again that 
\/2n z 2 +5e _2 e J ' z ' satisfies the fi = 2 root identity without need of any further 
correction factor e az . Since this is true for arbitrary z$ and without need of 
renormalisation, so it follows also for all fi G Z> 3 and thus Stirling's theorem is 
finally proven. 

One comment is worth noting regarding this constructive argument. Stir- 
ling's theorem actually expresses the asymptotic behaviour of T(z + 1) only as 
z — > +oo. The argument above, however, is naturally "two-sided" and suggests 
that r(z + 1) is asymptotically given by \/2~k z z+2 e~ z e J ^ also as z — > — oo, 
since in picking out the form of functions in £ so as to match the root side 
of our identities on the derivative side, we have had to call on functions (like 
z In z — z, In z, - etc) which are defined on all of R, not just (0, oo) or (— oo, 0). 

To check this, we can use the logarithmic version of the functional equation 
for r in (17) to deduce the asymptotic behaviour of ln(r(z + l)) as z — > — oo from 
its known behaviour as z — > oo from Stirling's theorem. We have, as z — > — oo, 

ln(r(z + l)) = (In7r + i7r) -ln(sin(7rz)) -ln(r(-z)) 
and, by Stirling's theorem, 
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-in(r(-z)) = -in(r(N)) = -Mr((N-i) + i)) 

>|-i)ln(|z|-l)-(M-l) + i^ 

1 1 I 

360 (|^| — 1) 3 + • • ■ 



+ 



1 1 



12 (z+1) 360 (z+1) 3 

f ( Z +|){(ln Z - W ) + I-il + Il- 



•} 1 



1 1 

' 360 z 3 



-(z+-)-in + 



+ -] 8 — ■ 

(z+i)lnz-z+ii 



l j_ 

360 z 3 



after simplification. Thus, overall, as z — > — oo, we have 



ln(r(z + l)) 



(z + \) In 2 — z + 



12 z 



1 1 

"360 z 3 



(Z + -) • i7T 



+ (lii7r + in) — ln(sin(7rz)) 



(43) 



We see that, as predicted, we do have the same asymptotic behaviour as per 
Stirling's theorem also when z — > — oo, except that in addition we have the 
linear and constant terms shown ((z + |) • in and (hi7r + in) respectively) and 
the term ln(sin(7rz)) which contains a countable collection of ln-divergences at 
the negative integer points. 

The reason our above derivation of Stirling's theorem nonetheless worked 
in spite of these extra one-sided terms as z — > — oo is that, owing to the 
factors on the derivative side of our root identities, the contributions made by 
these terms on this derivative side for /i € Z<o are 0. We shall omit calculations 
demonstrating this here, but note that in general terms this reflects the fact 
that it is asymptotic behaviour within a log-scale, rather than local divergences 
(even a countable collection of them), that matter in the computation of the 
derivative side of the root identities for \x e Z< . 
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4 The Generalised Root Identities for ( 



In the case of £, of course, we do not know the location of the non-trivial 
roots in advance and so our perspective must change. Unlike for T, where we 
calculated values for the root-sides of the generalised root identities by Cesaro 
means and used these to infer information about asymptotic behaviour of T 
from the derivative sides (e.g. Stirling's theorem), we now aim to calculate 
the derivative sides of the generalised root identities for /j G Z<o and thereby 
to investigate the location of the roots of ( by utilising Cesaro methods on 
the corresponding root sides, noting in particular that these Cesaro methods 
naturally depend critically on the geometric location of summands. 

4.1 The Derivative Sides of the Root Identities for ( 

In this case, for (, it turns out that there is no need to do detailed calculations 
on the derivative side using the Fourier definition; all that is required is the 
property (28) of (g^)^- Recall the Euler product formula for (, namely 

c(s)= n a-p" 8 ) -1 (44) 

pprime 

which is convergent for Jt(s) > 1. It follows that 
MCW) = - E ln(l-P~ S ) 

pprime 

= E {/' ~ • \p " • ',/' 35 • •••} (45) 

pprime 

and so, in light of property (28), we have that for 5i(so) > 1 the derivative side 
of the root identity for Q at \l is given by 

( p- So +2^- 1 p - 2s « +3"- 1 jr 3fl0 ) 

d c ( So , M ) = - Yl (^PTl \ (46) 

1 W pprtme [ +4^" 1 p- 4sa +... J 

This expression is clearly convergent for arbitrary /i for jft(so) > 1 an d so gives 
an expression for the derivative sides of the root identities for £ for arbitrary 
/jeC and jR(s ) > 1, which can then be extended also to Sft(so) < 1 by unique 
analytic continuation. 

Since, for arbitrary 3ft (s ) > 1 the sum in (46) converges to some finite value 
as \i — >• n for any n e Z< , and since T has simple poles at all the non-positive 
integers, so a corollary of the result in (46) is the following: 
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Lemma 4: When /i £ %<,q the derivative sides of the generalised root identities 
for £ are all identically zero as functions of So; that is 

rf(s ,M)=0 Vs V/i = 0,-1,-2,... (47) 

Since £ satisfies the generalised root identities, we should therefore have also 
r c (s ,M)=0 Vs Vp = 0,-1,-2,... (48) 

4.2 The Root Sides of the Root Identities for £ for fx = 

0, -1 and -2 

We now try to calculate the root sides of the generalised root identities for ( for 
the cases [i = 0, — 1 and —2 in order to verify (48) in these cases. 

4.2.1 Case (a): y, = : 

First we consider the trivial roots T. For so = we get the Cesaro sum of l's 
placed at the points 2,4,6,. . .. By dilation invariance, this equals the Cesaro 
sum of l's placed at the points 1, 2, 3,. . . and this is familiar from [1] as giving 
((0) = — |. For s = 2k this becomes instead the Cesaro sum of l's placed at 
the points 2k + 2, 2k + 4, 2k + 6,. . . which thus immediately leads to the value 
— I — k, and so in general 7 

]T Mifo, -1-0° = -^-^0 (49) 

Next, from the simple pole at s = 1 we get contribution —1 (since Mj = — 1 
here) . 

Now, for the contribution from the non-trivial roots, NT, in the critical strip 
we recall the explicit form of the Riemann-von Mangoldt counting function, 
N(T), which counts roots pi = fli + i-fi with < -fi < T. By Karatsuba-Korolev 
[5], this has the explicit form 



N(T) = N(T) + S(T) + -S(T) (50) 

7T 



where 



T T T 7 

^-s'-'s'-s + s (51) 

and S(T) is the famous argument of the zeta function (see e.g. [5] for formal 
definition), and 



7 The reasoning in going from so = 2k to arbitrary so to derive (49) here is deliberately 
heuristic, for reasons of brevity, but the result & similar deductions for the trivial root con- 
tributions when fi = — 1 and /i = — 2 are easy to deduce rigorously by Cesaro means if desired 
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7 ( ^^tan^m^r ^- {U} 1 



8(T) = - In 1 + — ~ + - tan" 1 — / ^ — L du (52) 

lJ 4 V 4 ^277 2i («+ 1)2 + (|)2 ^ > 

with {u} := u — Floor (u) being the saw-tooth function which rises linearly from 
to 1 on each integer interval [k, k + 1). 

Clearly for the p = root identity we are interested in finding the Cesaro 
limit of N(T) as T — > 00, at least after combining with the corresponding 
calculation for the roots below the real axis; i.e., writing T for the parameter 
tracing our count of roots with negative imaginary parts, we have 

Mi{s - Pif = Clim{N(T) + N(f)} (53) 

{so-NT} 

Now, by [5] we know that 5(T) = O(^), and by either [5] or Gordon-Sabeh in 
[6] we also know that if we define <Si(T) := J Q T S(t)dt in the usual way, then 

S 1 (T) = 0(\nT) (54) 
so that P[S]{T) = -> and ClimS(T) = 0. Thus (53) reduces to 

Mi{s - Pl )° = Clim{N{T) + N(f)} (55) 

{so-NT} 

As we have seen, however, in considering this Cesaro limit, we need to be adding 
in 1 for each root counted, not at the root pi itself but at s — Pi , and this 
geometry will be critical in our computations, as will the need to combine N(T) 
and N(T) before taking Cesaro limits (in particular we will shortly need this to 
resolve an apparent paradox between results for p = and p = — 1). To reflect 
this geometry we thus write 



so that 



s -(^+iT) and z = s - - if ) (56) 



T = i(z - (s - ^)) and f = —i(z — (sq - ^)) (57) 



and therefore 



and 



lnT^^ + ln,- ( ^-i ( ^- l( ^-... (58) 
2 z 2 z 2 3 z 3 y ' 



2 z 2 z 2 3 S 3 y ' 
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Here we are for the moment assuming that, at least for purposes of calculation 
of Cesaro limits, we may treat the parameters z and z as running purely along 
the critical line in enumerating the non-trivial zeros. We shall discuss this 
assumption (which would of course be justified if the Riemann hypothesis is 
true, and may be justified in any case by symmetry of roots around the critical 
line) further later. 
Thus we have 



T T 

— ln( — ) 



In — 4- i- — 



(so-A) l( S o-i) 2 



2 z 2 



'3 z 3 



'2n' 2tt 2 x 2n' 2tt 



v 2tt 



and therefore 



(so-h) .(sa-h) 



2tt 



■o(l) 



— M — )- — + ? = «( — )M — )-(-+*)( — )-— ( g o- -)ln( — ) 
2tt V 27r y 2f 8 V 2tt V 2tt ; V 2 A 27r y 2tt V0 2 ; V 2tt 

+(^o + ^) + o(l) (60) 

We see that taking care to place the NT roots relative to so by considering z 
rather than T introduces the extra terms of jSq and — | (which cancels the | 
to | ) in the constant term. If we add an equivalent constant-term contribution 
from N(T), arising from the roots in NT below the real axis, we get a constant 
term of ^sq + |, precisely as required to cancel the contributions from the trivial 
roots in (49) and from the simple pole and leave the root side of the /x = root 
identity with value as hoped - i.e. 

rc(*o,0)= ]T ^(50-^ = 

{sq— roots of 

in agreement with (48). 

However, there is a subtlety here and we must be careful. If we consider (60) 
alone then, recalling the dilation invariance of Cesaro limits (so that z -2- 
leaves limits unchanged) and the fact that z\nz and z are both pure Cesaro 
eigenfunctions (with eigenvalue |) and thus have generalised Cesaro limit 0, it 
follows that 
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But the term — ^(so — l) m (^) nas a P ure log-divergence which means 
that, on its own, it has no generalised Cesaro limit. Here, however, we see again 
the importance of being careful in both the geometric placement of roots and in 
considering the sum of N(T) and N(T) before taking the Cesaro limit, with T 
critically being treated as independent of T. Specifically, if we carefully mimic 
the calculations for (60) for N(T) in terms of z instead, taking account of the 
sign changes in (57) and (59) vis-a-vis the corresponding results for z, we end 
up with the corresponding relation 



L\ n (L^ L + l = _ i (_L) ]n (_L)_(![_ i )(_L) + _L( ao _I)i n (_L) 

2ir y 2Tr' 2tt 8 K 2Tr J V 2tt ; V 2 a 2tt ; 2tt V 2 7 V 2tt ; 

+(\so + l) + o(l) (62) 
In (55) we thus have, on noting (61) and its analogue for z, that 



E ^-^^r^l-Sizii^-lHd] (63) 

{s„-NT} ' K 

As discussed, on combining with the results for T and for the simple pole, we 
will thus have the desired result that, for fi = 0, 

r c (s ,0)= Mi(a -n) =0 Vs 

{so— roots of 

if and only if we have that 

Clim In ( |) = (64) 

z.z— 7-oo \Z/ 

But this is true in spite of the non-Cesaro-convergence of either In z or In z in 
the 1-d Cesaro setting because, as a 2-d function, ln(f ) can be ascribed a Cesaro 
limit and this limit is 0. To see this we first observe that f(z,z) = hi(|) is a 
direct eigenfunction of the 2-d Cesaro operator, P 2c i, with eigenvalue 1, since 



P 2d [ln (I)] = Pi d [lnz] - Pi d pnz] = (lnz - 1) - (In 5 - 1) = In (|) 

As such it can have a 2-d Cesaro limit, although the value of this limit is not 
directly clear. But then further note that since hi(|) and ln(-) should have the 
same limit, by symmetry, yet ln(|) = — ln(|), so this limit must in fact be 0. 

It follows at last that we have verified (48) and the generalised root identity 
for £ in the case n = by direct computation on the root side from the known 
form of N(T) in (56) - (59). This confirms that in this case the Cesaro count of 
all roots (T, pole and NT) of ( is identically for all s and has been achieved 
by showing that 
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]T Miiao-pi? = \*o + \ (65) 

{s„-NT} 

in a Cesaro sense. 8 

In passing we acknowledge that the above argument regarding (64) is some- 
what heuristic. This is discussed further in section 4.3 and [8], but we do note 
that, if desired, we could change perspective and instead view (64) as an impli- 
cation of £ satisfying the fx = root identity (rather than seeking to use (64) 
to verify this). This would at least then allow us to legitimately use (64) again 
in our next calculations for the \i — — 1 and fi = — 2 root identities, where it 
reappears. 

We thus turn now to the case of /i = —1. We will see that a naive calculation 
on the root side leads to an apparent paradox unless we are careful in again 
taking account of the geometric locations of z and z and of their independence. 
But doing so we can simultaneously resolve this paradox and verify that, on the 
root-side of the p, = — 1 root identity for ( we also have 

r c (a ,-l)= Miiso-n) 1 =0 Vs (66) 

{so — roots of 

as claimed in (48), at least modulo an estimate on S(T) which certainly holds 
conditional on RH. 

4.2.2 Case (b): u = -1 : 

In this case, for the trivial roots T, for s — we get the Cesaro sum of 2, 4, 6, . . . 
placed at the points 2,4,6,.... By dilation invariance, this equals the Cesaro 
sum of 2, 4, 6, . . . placed at the points 1, 2, 3,. . . and this is familiar from [1] as 
giving 2((— 1) = — g. For s = 2k this becomes instead the Cesaro sum of 
2k + 2, 2k + 4, 2k + 6, . . . placed at the points 2k + 2,2k + 4, 2k + 6,. . . which 
thus immediately leads to the value — | — k 2 — k, and so in general 

£ M^so-n)^-^- 1 ^- 1 - (67) 

{so-T} 

From the simple pole at s — 1 we get contribution 

M^so-n) 1 =-l-(*o-l) = l-*o (68) 

{so-pole} 

Thus, overall, from the combination of the trivial roots and the simple pole 
we get a total contribution to the root side of the ^ = — 1 root identity for ( of 

£ M I ( So -ny = - 1 -sl- 3 2 s + 5 - (69) 

{so-TUpole} 

8 Note as an aside that it follows trivially from the fact that the count of non-trivial roots 
in (65) is generically non-integral that there must be infinitely many non-trivial roots of £ 
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Turning now to the contribution from the non-trivial roots we have 



53 Miiso-Pi) 1 = J2 Miiso-fa-hi) 

{s -NT} {sa-NT} 

53 Mi(s -Pi)-i- 53 M ^ 



{so-JVT} {so-JVT} 

2' 



= {so-h 53 M i-i- 53 M ^ 



{so-NT} {so-JVT} 

5 3 

2~ u ' 4 So_ 4 



1 4 + 1^0 -I) -i- ]T M tll (70) 

{s„-NT} 



where here we have invoked the symmetry of NT roots w.r.t. the critical line to 
take (sq— Pi) outside the summation as (s a — |) and then utilised result (65) from 
the previous = computation. In these terms, the paradox alluded to earlier 
amounts to claiming naively that, because NT roots also occur in conjugate pairs 
symmetrically w.r.t. the real axis, so for real sq their contributions should cancel 
the remaining sum on the RHS in (70) to 0, thus leaving J2s So _ NT \ -^«( s o — 
Pi) 1 — (| s o + | s o — f )• This does not cancel the contribution in (69) to Vso 
and would thus leads to a contradiction of (48) and the generalised root identity 
for £ when \i = — 1. 

However this paradox is resolved by recalling that z and z must be treated as 
independent; and since the contributions to the sum in (70) from roots above the 
real axis and from roots below it are both individually classically divergent, we 
must work very carefully within a Cesaro framework to deduce the true Cesaro 
sum of X^{s -iVT} ^f»7»! rather than rely on naive pairwise cancellations. When 
we do this we will find that in fact the contribution from J2{ So -nt} ^*7» ^ s 
exactly as required to validate that, overall, on the root side of the fj, = —1 root 
identity for £ we do have 

r c (a ,-l)= 53 M i (so-r i ) 1 =0 Vs 

{so-roots of (} 

and thus (48) and the generalised root identities remain true for /i = —1. 

To see this, first note that for roots above the real axis we have the partial 
sum function of J^M^ given by udN(u) and similarly for roots below the 
axis. Thus, after accounting for the difference in sign of 7$ in the two cases, we 
get in (70) that 



53 M i {s -p i ) 1 = Q s o+ ^0- ^-iCZim udN(u) - ud7V(u)j 

(71) 



{so-NT} 

Now we shall adopt the notation that 
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S (T) := S(T), 5!(T)= [ T S (t)dt, S 2 (T) = f S 1 (t)dt , . . . 

Jo Jo 

and adapt it identically to define 

N (T) := N(T), and N,(T) = [ A^_i(t)di Vi G Z >0 , 

Jo 

N (T) := N(T), and N^T) = [ N^^dt Vi G Z >0 , and 

Jo 

<5o(T) := J(T), and S^T) = [ S^^dt Vi G Z >0 

JO 

for the functions introduced in (56) . 

In this notation, in (71) we get (on noting that all of N, N, S and 6 have 
finite limits as t -> 0+, with N(0) = 0) that 

f udN{u) = TN (T) - Ni (T) 
Jo 

= {TN (T) - N^T)} + {TS (T) - S^T)} + 1 {T5 (T) - ^(T)} 

(72) 

and similarly for udN(u). Now 

(i) Suppose we have the estimate that 

S 2 (T) = o(T) (73) 

Then since, in general, 

P[S n ]{T)= S ^± (74) 
it would follow immediately that P[S\\(T) — > as T — > oo so that 

CiimSi(T) = (75) 

Also, using integration by parts and noting Si(0) = by definition, we would 
have 

P[tS (t)](T) = i j T tSo{t)dt = S 1 (T) - S 1 (t)dt 
= Sx{T) - P[Si](T) 

4 (76) 
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Thus, combining (75) and (76) it follows in (72) that, conditional on the estimate 
(73), we have no contribution from the S'(T)-related terms, i.e. Clim(TSo(T) — 
Si(T)) = 0, and thus 

f udN(u) = {TN (T)-N 1 (T)} + -{TS (T)-S 1 (T)} (77) 
Jo 7r 

and similarly for J* Q T udN(u). 

Unfortunately, an unconditional estimate of the form (73) is not known for 
^(T). We do know, however, that a much stronger estimate for all S n (T) 
is in fact known, conditional on the Riemann hypothesis (RH), namely (see 
Titchmarsh in [7, pg 354]) 

Sn(T) = 0(^^) Vne^o (78) 

As such, we will proceed here on the assumption of (73) and (77) and thus 
end up obtaining a calculation of r^(s ,— 1) conditonal on (73), which is itself 
implied by RH. 

(ii) From the definition of N(T) :— ^ ln(^) — ^ + | we clearly have 

^>-{(£K>-(£) 2 +K£)} <™» 

and 



1 / T \ 2 , 7\ 3 / T V 7 / T \ 



2 " ! S '^'-i^ + SU <80) 



Thus 



™„(T)-tf l( T) = 2* - - m ( -)- 4 ) ) (81) 



and in light of (56)- (59) we get 
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f ^(^(*o-±)) 2 >^+*f-^-i^-...} 



2tt < 



2tt < 



2 \2ttj 



(«o-§: 



2?r \2ttJ \2ir 



: M^) + Hf + l)(^) 2 + ( 

+ (^)"f(s)-;( ! tf)Mi) 
+ (-i(W-'f-i( i ?#))+o(i) 



)M£) 



(82) 



In similar fashion, we have 



= 2tt < 



' -\m 2 Hi)+{^+\)m 2 +(^ 

-(^)-<!(£)-§(^)M£) 



2tt \2ttJ '"V27T 



+(+KW-<*-KW) +0 <*> 



(83) 



It follows, on taking Cesaro limits after combining TNo^-N^T) and TN (f)— 
Ni(T) in (71) (and recalling Climz 2 \nz = = Climzlnz = . . .), that we have 



Clim 

z,z—>-oo 



[TN (T) — Ni (T)] 



2.{(^#).Hf-,f)} - 

+Ste(-i( fi SJ e )'»H)) I 



(84) 



on again invoking the 2-d Cesaro argument used before to justify taking Clim In i 
as per (64). 
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(iii) This leaves just the terms from 6q(T) and S\(T) in (71) and (77) to 
resolve in order to complete the calculation of the root side of the ^ = — 1 root 
identity. 

Using the definition of S(T) in (52) we note from [5] that 



f 

Jo 



and in fact 

S(T) = ^ + ^ + ... (85) 

It follows that, in order to understand the Cesaro limit of TS(T) we need to 
first calculate a\, which is given by 

11 i ,. L 2 r (s-M) 



°i = 7^ + o - o T 2 77. \» T. d " (86) 



16 8 2t^x> 7 ( u+ i)2 + ( 



with the initial terms arising from the obvious expansions of In (l + j^j) and 
tan -1 (^f)- To calculate the last term in (86), although it is a classical limit, 
we will use a formal Cesaro computation and will likewise use such Cesaro 
methods in all calculations related to 5 in the remainder of this paper. Since 
these methods are formal and non-rigorous (see further discussion in section 4.3) 
we note immediately, however, that all results related to 6 are in fact rigorously 
verified, either theoretically or numerically, in Appendix 5.3. We start with the 
partial integral and apply a formal Taylor series expansion of the integrand: 



- fc +» (i - M) 



du 



4 

2^2 



4 

J"2 



— (^W du 



f k+a e,-{u})du 
^S» +a {\-{u}){u+\) 2 du 



rp2 



rnl a -l a2]+() 



i 

1 

ji4 



On recalling from [1] that Climot 11 = in (86) we therefore get that we have 



3 „ 1 1 

ai = 2 • — = — 

16 12 48 



(87) 



and thus, in (85), we have 
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<5(T) = ^ + |f + ... (88) 

This suffices to allow computation of the Cesaro limit of TSo(T), and would 
in fact suffice to allow computation of the contribution from (5-terms in (77) 
and (71) (since the constant of integration terms from Si(T) in the T and T 
contributions will cancel); but since we will need §2{T) in detail for the fj, = — 2 
root identity later, we shall calculate <$i(T) more carefully here anyway. 

From (52) we have 



Si(T) = 



In 1 



4t 2 



dt- 



2Jo 7o 



(«+i) 2 + (§) 2 



tan 



dudt 



2t 



dt 



(89) 



Now 



' In f 1 ■ )dl 



4t 2 



-f- 



+ 4t 2 



dt 



1. 



: T In 1 



4T 2 y 



1 

32 



In ( 1 + AT Z 



and 



2io 1 



ir. tan- m+ir^d* 



-T • tan" 1 — + — ln(l + 4T 2 ) 
In the final term, on reversing the order of integration we get 



11 (u + l)* + >i 



f 

Jo 



( 2 -W) 



H («+*)» + V 



dt dii 



du 



= < 



' /~(i-{ u })ln((«+i)2 + ^)du ' 
-2/ °°(i-{«})Iii(«+J)dti 
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Now the first of these two integrals can be expressed asymptotically in T as 

f(i- W ,{ 2 M|) + ^-...} d „ 

and a Cesaro argument identical to the one just used in deriving (87) (as one 
would expect) shows this is equal to 

I]n(T)-£ln2 + 0(-L) 
while the second of the integrals is a constant which we label A, i.e. 

A:= f\ l --{u})\n{u+\)du (90) 

and do not calculate further at this point. 9 

Thus, overall, combining terms we get finally in (89) that 



f lT 2 ln(l + 3 i 7 ) + ; |i n (i + 4r 2 ) + iT.tan- 1 (^) ) 
h(T) = { - /o °°(I- M) ln(( u+ I) 2 + f)d« 

+2/ oo (i-M)ln( U +i)d U 



where 



Combining (88) and (92) in (77) we then get 



(91) 
(92) 

(93) 



- {TS (T) — 5\(T)} 

7T 



and similarly 

l{fS (f)-S 1 (f)} 



■— In(T) + (— -CA+OiX) 
48 V ' V 48 ; ( T 2 ' 



1 

"48 



■ 7r 

\nz + i- 







■ - ,7T- 
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+ (^-Ci) + 0(^)} (94) 



+ (^-Ci)+0(J)}(95) 



9 We show in Appendix 5.3 that the integral defining A does indeed have a well-defined 
Cesaro value A ~ —0.104 
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It follows thus in (71) that 



Clim < 



^{T6 (T)-6 1 (T)} 
-ijTJom-JxCT)} ^ 



= CZim < — — — 



1 



48 48tt 



In 



2 

~48 



(96) 



on invoking the same 2-d Cesaro argument as before. 

Finally, combining (77). (84) and (96) in (71) we conclude that the root side 
of the (i = — 1 root identity for £ has contribution from NT of 



{so-NT} 

1 2 3 5 

= 4 So+ 2 S °-6 

and on combining with (69) this yields that, overall, 

r c (a ,-l)= ]T M^ao-rO^O Vs (98) 

{so — roots of 

That is, we have verified as promised that, at least conditional on estimate (73) 
for S2(T), which is itself implied by RH, the root side of the /j, = — 1 root identity 
for ( is indeed identically Vs as claimed in (48). 

Hence we now turn to the case of /i = — 2 and again try to perform a careful 
Cesaro calculation of the root side of the /j, = — 2 root identity for ( to see 
whether any insight can be gleaned from the claim in (48) that this, likewise, 
must be identically equal to 0. 

4.2.3 Case (c): y, = -2 : 

In this case, for the trivial roots T, for sq = we get the Cesaro sum of 
4, 16, 36, . . . placed at the points 2, 4, 6,. . .. By dilation invariance, this equals 
the Cesaro sum of 4, 16, 36, . . . placed at the points 1, 2, 3,. . . and this is familiar 
from [1] as giving 4£(— 2) = 0. For sq = 2k this becomes instead the Cesaro sum 
of (2k + 2) 2 , (2k + A) 2 , (2k + 6) 2 , . . . placed at the points 2fc + 2,2fc + 4,2fc + 6,... 
which thus immediately leads to the value — 4^^ =1 j 2 = — |fc 3 — 2fc 2 — 
and so in general 

J2 M i {8o-r i ) 2 = ~8l-±t%-±8 (99) 
From the simple pole at s — 1 we get contribution 



16 48' 

(97) 
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E Miiso-n) 2 = -1- (s -l) 2 = -s 2 + 2s -l (100) 

{s -pole} 

Thus, overall, from the combination of the trivial roots and the simple pole 
we get a total contribution to the root side of the [i = — 2 root identity for ( of 

E M i (s -r i ) 2 = -^s 3 -^s 2 +*s -l (1Q1) 

{s -TUpole} 

Turning now to the contribution from the non-trivial roots we have 
J2 M t (s - Pl ) 2 = E Mi{s - fa - i^) 2 

{s„-NT} {so-NT} 

E M,( So -ft) 2 -2z E M *( s o-ft)7* 

{s„-NT} {s„-NT} 

- E ( 102 ) 

{s„-NT} 

Writing = \ + ej the first of these sums becomes 

E m^-^f = E M t {so-\f+ E M * e * 2 

{so-JVT} {s Q -NT} ' {s„-NT} 

= ( S o-\) 2 (lso + l)+X £ (103) 

where 



X e := E Mid 2 = Climl^ e 2 (t)dN{t) + j e 2 (t)dN{t) 1 (104) 



on noting that the cross-term in vanishes because the non-trivial roots occur 
in mirror pairs either side of the critical line (and this holds identically for any 
given T so there are no Cesaro concerns with this cancellation). 

The sum in the second term becomes, on noting again that in going from T 
to T + dT any NT roots off the critical line occur in mirror pairs: 

E Mi(*o-A)7i = ( s o-^ E M ^ 

{s -NT} {s -NT} 

= {so-l)-i-(-\4 + \ao-^) (105) 

on using our results just computed for the non-trivial roots from the case of the 
fi = — 1 root identity (specifically (84) and (96)). 
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And the sum in the third term becomes 

i-T 



Ma?=Climl^ t 2 AN{t)+ PdN{t) I (106) 



{so-NT} 

Thus, overall, in (70) we get, after simplification, 



^ M t {s - Pi ) 2 

{so-NT} 



Now we have 



(7 „2 43 „ , n\ 
- Clim ( f n T i 2 diV(t) + f n T PdTV(i)} + Jf e 



(107) 



t 2 dN{t) = 



T 2 N a (T)-2 [ tN (t)dt 
Jo 

T 2 N (T) - 27Wi (T) + 2 f Ni (t) dt 

Jo 

T 2 N (T) - 2TN l (T) + 2N 2 (T) 



on noting iVi(O) = N 2 (0) = q ; i.e. 



/ t 2 dN(t) = T 2 N (T) - 2TN X (T) + 27V 2 (T) 



(108) 



and similarly for J Q T t 2 dN{t). 

Recalling the original expression for N(T) = N(T) + S(T) + ^S(T) it follows 
in (107) that we have 
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Clira ■ 

z.z—xx, 



= Clim < 

z.z—xx, 



t 2 dN{t) + J^ PdJV(t)| 
|T 2 iV (r) - 2TN 1 (T) + 2N 2 (T)} 
+ [f 2 N (f) - 2XWi(T) + 2iV 2 (T)} 
+ |T 2 5 (T) - 2T5i(T) + 2S 2 (T)} 
+ [f 2 S Q (f) - 2TS 1 (f) + 2S 2 (f )} 
+ i {T 2 S (T) - 2T6 1 (T) + 2S 2 (T)} 
+ i {f 2 S (f) - 2T6 1 (f) + 25 2 (f)} 



(109) 



Now 

(i) As in the case of [i 



-1, suppose here that we have an estimate 



S 3 (T) = o(T) (110) 

which, as before, would certainly follow from the much stronger estimates (78) 
if RH is true. Then we would have that ^ J Q T S 2 (t)dt = ^iZ) = (i) so that 
P[S 2 ](T) and thus ClimS 2 (T) = ClimS 2 {f) = 0; and, given the capacity 

z— >oo z— >oo 

to turn Cesaro limits of T 2 So and T5i into Cesaro limits of S 2 (in the same 
fashion as per the argument in (76)), it would follow once again that there was 
no contribution in (109) from the S'-terms. 

As in the case of /i = — 1, no such estimate (110) is in fact known uncondi- 
tionally, but we shall proceed here on the assumption of (110), so that 



The contribution of the S — terms to (109) equals (HI) 

and hence continue on towards a calculation of r^(s , —2) conditional on (110), 
which is itself implied by RH. 

(ii) From our earlier computation of Ni(T) in (80), on letting u — ^r, we 
have 
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and therefore, on combining (57), (80) and (112) and simplifying, we get that 



T 2 N (T) - 2TN 1 (T) + 2N 2 {T) 



Thus 



T 2 N (T) - 2TN 1 (T) + 2N 2 {T) 



(2tt) 2 



{z-{s a -\)Y 

(27T)3 



In 



Z I ' 7T 



2 z 2 3 z 3 

. 1 1 



+ 



9 ^ (2tt)3 ) 



and similarly 



T 2 N {T) - 2TN 1 {T) + 2N 2 (T) 



= (27T) 2 



3 ^ (2tt)3 ) 



i («o-|) 2 _ i («o-|) 3 _ 

2 z 2 3 z 3 

'(z-( So -i)) 3 ^ 



i f(z-(s -h)) 3 \ 

9 { (2^)3 J 



(113) 



(114) 



(115) 



Thus, bearing in mind our usual observations that Clim z In z = etc, we have, 

z— >oo 

after omission of these terms and simplification, that 



Clim 

z.z— >oo 



Clim 

z.z— >oo 



|T 2 iV (T) - 2TN 1 (T) + 2N 2 {T)} + ' 
|f 2 iV (f ) - 2TW 1 (T) + 2/V 2 (f )| 
-^K-|) 3 + 



-6^-2) 3 



IK -i) 3 -^ In I 
1 



(116) 
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on invoking the same 2-d Cesaro arguments as before, 
(iii) From (91) we have 



' §/ Vln(l + ^)d;+JJ>(l + 4t 2 )dt - 
5 2 (T) = l +i/ ? Ytan- 1 (i)dt 

, -J T / OO (|-M)ln(( W +|) 2 + $)d U di + 2AT 



(117) 



where A is defined as per (90). Now 
(a) 



t 2 ln 1 



At 2 



dt 



(b) 



1 T 3 
8 ' ~3 

1 T 3 
8 ' T 



In 



1 + AT 2 ) + nJ Q i + 



1 + At 2 



at 



In 1 



1 



1 



1 



AT 2 7 48 96 



T--tan" 1 (2T) 



(118) 



32 



/ 



In (1 + At 2 ) At 



(c) 



32 
3 

32 



T In (l + 4T : 



o 1 + ^ 2 



dt 



T In (1 + AT 2 ) - At + ^ tan- 1 (2T) 



(119) 



i / t -tan" 1 ( -J- ) dt = l-T 2 tan- 1 ( J=) + 4 / 
4 7 ^2*/ 8 \ 2T J 16 l + 4t 2 



4t 2 



^tan- 1 



1 



dt 



(120) 



and 



(d) Reversing the order of integration we get 

t 



/•°° 1 



4 ) 2 + -)dM. 



t 2+4(«+i)2 



dt ^ du 
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,oo 1 

= L «5- W) 

/•°° 1 



' Tin ((„+!)=+?) -2T 

k + 8 ( u + i) 2 /o T ranF di 

' Tln(( M +i) 2 + ^)-2T 

+4( U +i)tan- 1 ( 5 ^ iy ) 
2(lnT-ln2) + i^^ + 



• du 



■ du 



- 2T 



K _ 2(u+\) 
2 T 



+4( U +i) 
2{/ ° (|-{«})d«}TInT 
(2 + 21n2){/ oo (i-M)d M }T 
k +2 7 r{/ oo (i-M).( U +i)d U }+O(i) , 



> du 



(121) 



Now 



(d)(i) We saw before that Clim J" — {u})du = ^ so that, in Cesaro 
terms, 



r°° i i 
/ <2" 12 



(122) 



and 



(d)(ii) Similarly 



pk+a 

L 1 



(x-{«})-(«+ 7 )d. 



fe-1 „i 



= E / 



1 f a 1 1 

a)(j + a + -)da + J (- - a)(k + a + -)da 



1 1 1 

~ 2 )da + (k + -) I (--a)da+ (-a - a 2 )da 
' " Jo 2 



4 / ./o v 2 



1 , / j lwl 1 2, 1 2 1 3 

12 fc+ ( fc+ 4 )( 2 a -2 a ) + 4 a 



= — -fca 2 + -fca - 



1 



1 



1 



1 



k + (--of + -ct + -a) 



(123) 



2 2 12 v 3 

But, extending lemma 12 in [1], we have in fact (see proof in Appendix 5.4) that 



Lemma 6: // f(k + a) = k n a r then 
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CJi S fik + a) = {-ir-± Tl (124) 
Applying this in (123) we thus get, after simplification, 

f k+a 1 1 1 

( 2 - {u}) - (M+ 4 )dU= 48 

i.e. in Cesaro terms 

/ °(^-M)-(«+J)d«=^ (125) 
Then, combining (122) and (125) in (121) it follows that we have, up to O(^), 

/•T /-oo -i / 1 j.2\ i i 

^ / o (2-W)ln^+ J ) 2 + T Jd U dt=-TlnT--(l + ln2)T+- (126) 

and therefore, finally, combining (118), (119), (120) and (126) in (117) we get, 
after simplification, 

5 2 (T) = ^TlnT + BT-^+Oi^) (127) 

where 

B = — + —ln2 + 2^ (128) 

96 48 y ' 

Combining (127) with our earlier expressions for 6q{T) and 8\(T) in 
and (92) it follows that we have, in the (5-terms in (109), that 



1 f T 2 6 (T)-2T6 1 (T) ) ( -LT-±TlnT-2C\T ) 

+26 2 (T) J n { +±T\nT + 2BT - § J J 

l{ DT -h) + °^ ^ 



where 



D= — -2d+2B = - — 
48 48 



and 



similarly for \ {t%(T) - 2T5 l {f) + 2£ 2 (T')}. 



On taking C7im we thus get contribution from these <5-terms in (109) of 
± g Kj» • (* - (-o - ^)) - Di ■ (z- (s Q - I)) - ^} = -1 (130) 
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on noting Clim z = Clim z = in the usual way. 

z— >oo z— >oo 

Finally, combining our results from (111) (conjectural or conditional), (116) 
and (130) in (109), we obtain in (107) that we have 

„ f (I s o ^ H s o + H) 1 

J2 M t (s - Pl ) 2 = I \ 

{sa-NT} [ +i( So -|) 3 + ^+^ e J 

= ^ S + ^0-^0 + ^+^e (131) 

And then, finally, combining (131) in turn with (101) we find that, at least 
conditional on estimate (110) for S^(T), which is itself implied by the RH, the 
root side of the [i = — 2 root identity for £ is given by 

r f (s ,-2) = M^-nf = - l -+X t Vs (132) 

{s -roots of C} 

From this and the fact that, according to (48) this should equal iden- 
tically, we can make the following claim (modulo the rigour of the preceding 
computations) : 

Result 1: The Riemann hypothesis (RH) is false - that is, there exist non- 
trivial zeros of £ off the critical line jft(s) = \ in the critical strip. 

Proof: We argue by contradiction. If RH were true, then 

(i) X t = trivially 

(ii) The conditional results regarding Cesaro limits of S 2 (T), Si(T), and 
So(T) noted in (111) all actually do hold, by the explicit estimates, conditional 
on RH, given in [7, pg 354], and 

(iii) All the NT roots pi are of the form exactly \ + iji (7$ € R) and so the 
use of z and z exactly to parametrise the counting functions N(T) and N(T) in 
C becomes legitimate, since these really do then become step functions on the 
critical line. 

But then, in (132) we obtain that 

r c (s ,-2)= Mi{s - nf = Vs 

{s -roots of (} 

which contradicts the requirement, from the fact that £ satisfies the generalised 
root identities (and in particular the fj, = — 2 root identity), that r^(so,— 2) 
should be identically zero for all s as per (48). 
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4.3 Discussion of Issues Requiring Further Attention 

The emphasis in this paper has been on the development of certain methods 
pertaining to remainder Cesaro summation/convergence, and on calculation us- 
ing these methods. Consequently, there has been a conscious de-emphasis on 
formal rigour and detailed proof. In light of this, and especially of the claim in 
result 1 that such computations imply the Riemann hypothesis is false, we now 
conclude by trying to clarify which claims are in fact rigorously supported and 
where there are gaps requiring more detailed justification. 

(1) The first major issue requiring more rigorous justification is that, in the 
case of (, we have simply asserted that because £ satisfies the generalised root 
identity for /z = 1 for all so (modulo the ir~i obstruction factor), so it must (like 

T) continue to satisfy these identities for all sq for fi = 0, —1, —2, In part 

this reflects the belief that, having verified the /i = 1 identity, so that roughly 
-f^y = T,{s -rootsofC} W=u) for arbitrary s o (the more rigorous statement 
including the renormalisation correction and 7r _ 3 term is (26)), so further dif- 
ferentiation by ^g^^ should imply that the generalised root identities hold 
for £ also at arbitrary /i, and in particular /i = 0, —1, —2, . . .. But as well as be- 
ing heuristic this skates over the possibility that there exist functions which act 
as obstructions for the root identities for /i = 0, —1, —2, ... in the same way as 
the functions e az act as obstructions for the root identities for [i € Z>i. Since 
n~i((s) does already satisfy the /j = 1 identity for arbitrary s (and hence for 
H = 2,3, . . .), so this is the question of whether, for any given n G Z< , there 
exists a function, g(s) which fails to satisfy the root identities for \i = n but 
satisfies them all at arbitrary s for fj, = n + 1, . . . 1, 2, 3, . . . (since we need to 
leave these identities undisturbed). We have not addressed this question at all 
in this paper although, based on considerable experimentation, we believe that 
the non-existence of such additional obstructions is a reasonable conjecture (as 
well as one that seems quite distinct from RH). 

Another way of thinking of this issue is that we know from [1] (e.g. the 
discussion of removable singularities for discrete Cesaro schemes discussed in [1], 
section 4.3) that the values obtained by Cesaro methods at individual parameter 
values (e.g. /i = 0,-1 and — 2) are only reliable if they arise from Cesaro 
analysis in an open region around these values which represents a proper analytic 
continuation (e.g. in /x) to these regions. Since we have only actually verified 
the root identities g^(so,m) = r^(so,n) f°r C at the isolated values /i e Z>i, 
we have not in fact even established a starting open region in C where the root 
identities are true, let alone systematically extended to regions encompassing the 
critical points /i = 0, —1 and — 2. In a separate paper ([8]), however, we address 
this deficiency by showing numerically that it appears that the generalised root 
identities for £ are indeed satisfied for all real /x > 1, where both d^(so, n) and 
r c( s 0i/ i ) are given by classically convergent expressions; and furthermore that 
they appear to continue to be satisfied both for values < fi < 1 where Cesaro 
divergences need to be removed on the root side and for values — 1 < \i < 
where Cesaro divergences not only need to be removed but a further Cesaro 
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averaging of the resulting residual partial-sum function applied. For /x further 
to the left (i.e. /x < — 1) the numerical difficulties in verifying the root identities 
then become too significant without major effort beyond that attempted in [8] 
owing to the need to perform additional applications of the Cesaro averaging 
operator P, but nonetheless these numerical results in [8] seem to us to be 
strongly suggestive that ( does indeed satisfy the generalised root identities for 
all /x G R (and hence all /x G C by analytic continuation) as desired. 

In addition to this we would also argue that the fact that we have computed 
the root sides explicitly when fj, = and /i = — 1 and thereby demonstrated that 
£ does satisfy the generalised root identities at these points (modulo estimate 
(73) in the case of it = —1) is strongly suggestive on its own that £ should 
continue to satisfy these identities at /i = —2, which is the key to result 1; 
and we believe this especially because it seems to us that the nature of the 
calculation in the /i = and fi = — 1 cases is non-obvious and itself suggestive 
of the aptness of the Cesaro approach. 

(2) A second major issue is that we have not, in fact, rigorously validated the 
definition of generalised Cesaro convergence along a contour given in definition 
1, i.e. we have not rigorously shown that this definition, and in particular the 
distinction between removing "eigenfunctions" in the geometric variable z = 7(f) 
but then just averaging in arc-length t, guarantees correct analytic continuation 
in the situations in which we have applied it (e.g. in z and fi on the root- 
side, r(zo,/i), of root identities). The numerous examples in the paper where 
this is confirmed by calculation (e.g. the definition of V, the root identities for 
T at it e Z< and for ( at fi = and fi = — 1) give strong "experimental" 
evidence that this is so, but no rigorous argument has been given. Briefly, such 
an argument would likely proceed along the following lines, illustrated here in 
the context of the simpler case of the root side of the root identities for T(z + 1) 
(but easily adaptable to the case of £) : 

(a) Initially restrict to zo real so that for rr (z n , /i) we are back in the setting of 
partial sum functions on [0, 00) handled rigorously in [1] via regular polynomials, 
q(P), in P. When ji is real and /i > 1 the series defining rr(zQ,fi) is classically 
convergent, and when /x < 1 the "removable" Cesaro eigenfunctions are all in 
fact naturally functions (e.g. z p etc) in the geometric variable z = z + k + a 

(b) Thus the rigorous methods of [1] should establish, for all zq real, the 
correct analytic continuation of rr(zo,M) from \i > 1 to all real /x. Analytic 
extension of r r (z ,M) to z off the real line (holding /x real still initially) is 
then uniquely determined and seemingly must involve still "removing" the same 
divergent terms in the variable z, which is no longer real but has the same 
geometric meaning, namely z = z + k + a as described in definition 1 

(c) Finally, analytic continuation to all /x off the real line is then also uniquely 
determined (although it is no longer necessarily still interpretable in terms of 
Cesaro convergence if z is not real; note that in this paper we have never 
actually considered non-real /x). 

With respect to (b) , note that the fact that we remove "eigenfunctions" such 
as z p in the geometric variable z (rather than, say t p in t = k + a) means that 
as ji varies and causes p to pass through positive integer values 1,2,... we do 
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not instantaneously pick up discontinuous contributions to the Cesaro limits at 
these [i- values ( as we would if we were removing "eigenfunctions" in t) . This 
observation bears directly on the remarks (ii) and (iii) made in commenting 
upon definition 1 in section 2.1, and is especially relevant for why we hope that 
we can still trust the Cesaro evaluations in this paper for V and ( at non-positive 
integer values without need for correction of removable singularities. 

Note finally, in passing, that the adaptation of the above scheme of argument 
to the case of ( would involve negligible change for handling the trivial roots, 
and would merely require an adaptation of the approach in [1] from functions 
on [0, oo) to functions on the positive and negative imaginary axes for handling 
the NT + and iVT_ components of the non-trivial roots. 

(3) For ( in the cases of /j, = 0,-1 and — 2 we have relied on a claim 
that we can ascribe a 2d Cesaro limit of to Clim In (f ). In fact, along with 

a parameter (say w) for the trivial roots, T, the calculations for r^(sQ,fi) at 
fj, = 0, — 1, —2, . . . are all really 3d Cesaro calculations, but we can treat the 
trivial root calculations independently and it is only in relation to the claim 
that Clim In (I) =0 that we are obliged to consider the two variables z and 

z,z— >oo 

z simultaneously. The argument for this claim (in (64)) is, however, clearly 
non-rigorous (particularly since we have not formalised the 2d notion, on a pair 
of contours, being used). Nonetheless, we make the following two observations 
about this claim here: 

(a) First we reiterate that, as noted in section 4.2.1, for purposes of our 
claim re RH this 2d Cesaro limit could alternatively be taken as an implication 
of the /i = root identity and then applied without caveat to the fj, = — 1 and 
H = — 2 cases on that basis, and 

(b) Secondly, in [8], as part of explicitly considering the Cesaro treatment 
of the non-trivial root contributions to r^(s ,^) in a neighbourhood of 0, we 
provide strong numerical evidence to justify taking Clim showing In (4) = 

z,zby— >oo 

that r^(s ,^) — > as [i approaches from below. Moreover we give a precise 
conjecture as to how this ln(§) divergence may be arising naturally in the 
limit as fi — >• (essentially from the contributions from NT + and iVT_ each 
individually diverging as \i — > but in opposite directions so that the overall 
iVT-contribution remains finite) and likewise give strong numerical evidence 
that this conjecture is true. This would validate fully that taking Clim In (4) = 

z,z— >-oo 

at u — (and presumably similarly at = — 1 and fj, = —2) is the correct 
thing to do in order to ensure correct analytic continuation. 

(4) Lastly in order to focus on calculation, in numerous areas throughout 
the paper we have been somewhat loose, for example in: 

(a) Adopting a heuristic approach in applying Fourier theory to perform 
calculations on the derivative side in some examples of the generalised root 
identities, and indeed even in regards to the Fourier definition of (jj)^ on the 
derivative side of these identities in the first place, 

(b) Interpreting integrals, or , in a Cesaro way as required (effec- 
tively applying a universal Cesaro viewpoint throughout the paper), 
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(c) Utilising Taylor series expansions without always being careful regarding 
their domains of convergence, and 

(d) Using formal Cesaro methods to obtain the required parts of the asymp- 
totic expansions of S, Si and 5 2 in the calculations of the root sides of the root 
identities for £ at \i = 0, — 1 and — 2. 

For instance, in analysing 6, Si and S 2 as mentioned in (d), despite the fact 
that these are all classically well-defined functions and thus should be calculable 
without recourse to Cesaro methods, we instead began by applying a Taylor 
series expansion 

1 4( M +i) 2 I6(u+\r 

1 i(u+\y t 2 T 4 

inside an integral, f™Z^° (namely the integral term in the definition of S) , where 
for any fixed T this expansion is not classically valid for most of the domain of 
integration < u < oo. We nonetheless proceeded by working formally term by 
term in descending powers of T, using Cesaro methods (specifically lemma 6) to 
evaluate the resulting divergent integrals in u which then form the coefficients. 
And, for S, we thereby ended up with result (88), which can readily be extended 
by keeping the next lower order terms in T to give a more detailed expansion 
such as 

6{T) = 48 T + 5760 + 80640 T^ + °(t^) (133) 

and similar results for Si(T) and ^(T 1 ). 

Notwithstanding this "looseness," however, we believe that all the calcula- 
tional claims in this paper are solid. In particular, in relation to the listed areas 
(a) - (d) and how they bear on our calculations for £, note that: 

(i) On the derivative side of the root identities for £ at arbitrary fi we have 
only ever relied on property (28) of (i)'* and not needed to explicitly invoke 
Fourier and distributional arguments 10 , and 

(ii) As noted in section 4.2.2, all the claims regarding S, Si and <5 2 used in 
calculating the root sides of the root identities for ( at /; = 0,-1 and — 2 are 
in fact true and have been verified both numerically and by alternative rigorous 
derivation (see Appendix 5.3) 11 . 

Overall, despite the issues (l)-(4) flagged in this section (some of which are 
further addressed in [8]), we believe, both for the variety of reasons mentioned in 

10 We have in fact attempted a derivation of lemma 4 using Fourier and distributional 
methods. The arguments rely on identifying the asymptotic behaviour of as s — > oo 
and hence also as s — > — oo via the functional equation; to handle the T-style asymptotic 
behaviour that thereby arises as s — > — oo we adapt our arguments regarding Stirling's theorem 
in section 3.4 to this "one-sided" setting (i.e. Stirling asymptotics only as s -> — oo not in 
both directions). Details are omitted here however. 

11 The reason we have nonetheless adopted the formal Cesaro approach in analysing these 
functions in the paper itself is partly for reasons of brevity, partly in order to keep a Cesaro 
perspective throughout the paper and partly because we believe that the fact that these Cesaro 
calculations work successfully despite their formalism is in itself interesting 
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the discussion in this section and on the basis of the results derived throughout 
the paper in a range of examples, that the methods and results of this paper 
are correct. 

5 Appendices 

5.1 Dilation Invariance of Cesaro Convergence 

The inverse of the Cesaro operator P[f](x) := - f£ f(t) dt is easily seen to be 
given by 

P- 1 = x^- + l (134) 
ax 

(as discussed in [1]). On the other hand, the generator of the dilation group 
D r : x h-> rx , r > 0, is given by x-^, since for e very small we have 

f((l + e)x) = f(x) + f'(x)-xe + 0(e 2 ) 

so that 

D 1+£ = I + e ■ x-^- + 0(e 2 ) 

It follows at once that dilations commute with the application of P and so 
Cesaro limits of functions must be dilation-invariant. 

A similar, though more involved, argument shows that, in fact, Cesaro limits 
are also invariant under rescaling (S r : x i-> x r , r > 0), but this has not been 
utilised in this paper and so we omit details here. 

5.2 Figures 

The first two of the following diagrams have been referenced in the body of the 
paper. 

Figurel: Cesaro Dilation Invariance and the Duplication Formulae for T 
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Figure 2: Bi-directional Remainder Summation and the Functional Equation 

r 
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A 




X 

The final figure below depicts relationships among the concepts and results 
of this paper, and linkages with other well-known results. It is included in order 
to provide a survey of the avenues investigated in the paper, and to illustrate 
why they may provide a new perspective on certain issues, in particular the way 
they combine in leading to the claim made in Result 1. 

Figure3: Connections Among the Results of this Paper 
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5.3 Rigorous Verification of Formal Cesaro Results for 5, Si 
and 82 

In the course of evaluating the root side of the root identities for £ when \i = — 1 
and \x = — 2 we used the following results regarding do (= £), Si and £2 which 
we derived using formal Cesaro means: 



5(T) 



11 7 1 31 1 

48 T + 5760 + 80640 



1 

T 5 



(135) 



= i lnT+ 4 + S ln2+2A >- n^^-3^^ 

and 



1 

T 4 



(136) 



^^ rlnT+ ^ + S ln2+2A > T -6l + nk^ + ^^ 



1/1 

(137) 



where here A is the constant given by the Cesaro value of the integral 

A = J™(±-{u})]n(u+±)du (138) 
These implied the relationships 
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(139) 

and 



T 2 5 (T)-2TS 1 (T)+2S 2 (T) = --T- — + — - + — \+o ( X) (140) 

uv ; 1 ; v 1 48 32 5760 T 48384 T 3 \T 3 J v 7 

which constituted the 5-term contributions to the roots sides of the /j, = — 1 and 
/i = — 2 root identities for £• 

We now show that these results are all in fact valid despite their heuristic 
derivation in the body of the paper. We start with a proof of the asymptotic 
expansion for 5 in result (135). 

Proof of Asymptotic Expansion for 5: Recall that 



(^ + i) 2 + (i) 



Using that ln(l + t) = t — \ + y — . . ., we first get the following expansion, 
which converges classically if T > \: 

T ( 1 \ _ 1 1 11 1 1_ 

4 n ^ + AT 2 ) ~ 16 T 128 T 3 + 768 T 5 

Similarly, using tan -1 (t) = t— y + y — ...we get the following expansion which 
also converges classically if T > \: 

1 _ x ( 1 \ _ 1 1 11 11 
4 tan yxf) ~ 8T " MT 3 + 640T5 

Next, write the final term in the expression for S as 

gri^B^du (141) 



/ 

Jo 



where 

T 
2 



9t(u) 



(«+|) 2 + (|) 2 

and Bi(u) = {u} — |is the first "periodic Bernoulli function" as defined in section 
1.1. Using the fact that £^B n (u) — n_B„_i(u) for n > 2 (except at the integer 
points when n = 2 which is a set of measure 0), we can write the integral (141) 
as 
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1 f°° d ~ 1 - If 00 

2J 9t(u) —B 2 (u) du = ^9t(u) B 2 (u)\^ =0 - -J g' T {u) B 2 {u) du 

and since gr(u) — > as u — > oo and ^(O) = i?2 this is 



-f 9 ,(0)-i 



<7^(w) B 2 (u) du 



2 



1 f°° d - 

9t(0)--J g' T {u)—B s (u)&u 



Now i?3(0) = B3 = and it is easy to check that jy(u) -> as u 
integrating by parts again, the integral (141) becomes 



oo, so 



B 1 /*°° 

-^3T(0) + -y o 5 ^( U )B 3 ( U )d U 

= -f.M0) + 24 / o ^(«)^(„)d U 

which, since g'^(u) — >• as u — > oo and £?4(0) = £>4, becomes in turn 

g'T(u)B 4 (u) du 



-%■((>) -&#(o) : 



24' 



24 



Continuing in the same way, noting that g!^\u) — > as u — > oo and that 
Sfe(O) = Bk with i?2i+i = for all I > 1, we obtain the following expression for 
the integral (141) 



f 

Jo 



gr{u) B\{u) du = < 



-f^(0)-f^(0)-ff <#>(()) 

+i/;#( u )B r ( I1 )d I i 



(142) 



which could be extended further if desired. 
In this case we have 
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g' T (u 
9t(u 



(4)/ 

9t \ u 



(5) / 

9t \ u 

(6) / 

9t \ u 



((-+|) 2 + (f) 2 ) 2 3 
T(3(,+ j) 2 -(f) 2 ) _ 1 
((-+i) 2 + (|) 2 ) 3 3 
12T(-( M +i)3 + ( M +l)(^) = J_ 
((«+i) 2 + (|) 2 ) 4 ~ 5 

12T(5( M +i) 4 -10( M +i) 2 (j) 2 + (j) 4 ) = 1 
((«+|) 2 + (|) 2 ) 5 5 
1ST (-24( M + i) 5 + 80( M + |) 3 (f ) 2 - 24( M + | Xj ) 4 ) = 1 
((«+i) 2 + (|) 2 ) 6 7 

45T (56( M + |) 6 - 280( M + |) 4 (f ) 2 + lggfo + j) 2 (| ) 4 - 8(f) 6 ) = 1 
( (u+ I )2 + (f) 2 ) 7 ^ 



But it follows from the expression for g^\u) that in (142) we have 

^ g¥\u)B r (u)du = o(-) 
which we shall show by demonstrating that 



f 

Jo 



g^\u) B*[{u) du 



satisfies It — > as T — » oo. To see this, note from the expression for g^\ u ) 
that 

(("+l) 2 +(l) 2 ) 4 

and so 

/■OO 

It = hT(u)Bj(u)du 
Jo 

where < for a constant C u depending only on u. Moreover, 

T 5 \9 { t\u)\< 



((-+|) 2 + (?) 2 ) 

and differentiating the right hand side with respect to T, we see that, for fixed 
it, the right hand side is maximised when T = C(u+ j). This maximum /i(u) is 
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a constant multiple of j^ry- Hence J °° h(u)du < oo and since \}it(u) Bj{u) \ is 
at most a constant multiple of h(u), so we can apply the Dominated Convergence 
Theorem and see that It — > as T — » oo. Thus in (142) we have that the integral 
(141) is given by 



-f^(o)-t^(o)-f^(o) + o( i,) 

2T 16T(4T 2 - 3) 512T(16T 4 - 40T 2 + 5) 1 



3(4T 2 + 1) 45(4T 2 + 1) 3 315(4T 2 + 1) 5 V T 5 ' 

11 7 1 5 1 1 
~6T + 360T 3 ~ 2016 + °'t^ 

Adding in the contributions from j In (l + and |tan _1 (^) we obtain 
the asymptotic expansion 

5{T) 1 1 7 1 31 1 ( 1 



48 T 5760 T 3 80640 T 5 \T 5 

exactly as claimed in (135). 

Next we verify that the integral defining A in (138) does in fact have a well- 
defined Cesaro value A ~ —0.104 as claimed. To see this, consider the partial 
integral 



= < 



(g - W)M«+-)dM 
E J to 1 /o(5-5){lni + ^ 



+ -a){]n*+<*+i> 



i) i(a+J) 2 



+ . . .}dc 



i(«+|) 2 



A: 2 



+ . . .}da 



+(|a- \a 2 )\nk + C + 0{\) 

3=0 J 

= -\^-a + \)Xnk +{ C-^) + { \) 



where 7 ~ 0.577 is Euler's constant. But since ^{a 2 — a + |)da = it follows 
easily that P[(a 2 — a + |) lnfc](fc + a) — > as fc — » 00. Hence we have that the 
partial integral J" — {«}) ln(u+ \)du does have a generalised Cesaro limit 



as k — >• 00 and in fact we have 
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= Km 2/ (--a){In(j+a+-)-Inj-^ T ^}da-— 7 
fe ^°°^o^o 2 4 j 12 

(143) 

which we may calculate numerically to get the above approximate value. 

This leaves just the verification of the asymptotic expressions (136) and (137) 
for 5i and 82- Fully rigorous derivations of these have also been obtained along 
the lines of the derivation above for 8, after obtaining recurrence relations and 
certain preliminary estimates for the integrals 

B 2r -i(u) 

-au 



(("+|) 2 + (f) 2 ) m 

For the sake of brevity, however, we omit details here. Instead, as an alternative 
we note that the expressions in (135)-(137) for So, £>i and 82 can also all be readily 
verified numerically by writing the integral in the definition of 8 as a convergent 
sum over j of integrals from j to j + 1; since on each such unit interval the 
integral is easy to perform exactly (and likewise for the furrther integrations 
for 8 1 and 8 2 ), so the resulting expressions for <5o, 8\ and 82 are comprised 
entirely of elementary functions and convergent sums of the same, allowing 
thorough numerical verification of the claimed asymptotic expansions. Such 
sample numerical confirmations, implemented via VBA code, can be found in the 
XL2007 spreadsheet "delta_0_l_2_AsymptoticsFinal2.xlsm" made available 
with this paper. 

5.4 Proof of Lemma 6 

In order to perform our calculations related to 8 in our derivations of the values of 
r^(s , —1) and r^(so, —2) we invoked the result (124) given as lemma 6, namely 
that 

Climk n a r = (-l) n 1 (144) 

fe^co n + r + 1 

In these cases the specific claims required (that Clim a r = -j-=- , Clim ka 2 = — \ , 
Climka = — = , and Climk = — ^) all involve only n = or n = 1 and are 

k— >oo k— >oo 

readily checked independently by direct calculation. However, since we regard 
lemma 6 as an independently interesting result in the Cesaro framework, we 
include its general proof here now. 
Writing x = k + a as usual, let 

fn,r(x) ■■= k n a r and g n ,r(x) ■= x n a r 
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When n = note that fo, r (x) = 5o,r(x) and it is trivial that they both converge 
in a generalised Cesaro sense to ^jpx, i.e. 

Clim for{x) = Clim go.r(x) = — — — (145) 

x— >oo ' ' x— >oo ' r+1 

We first prove the following additional lemma: 
Lemma 7: 

c (-1)™ c 

fn,r( X ) -> — 7" if an d onl y if 9n,r{ x ) ~> ^ n <= Z >0 (146) 

?2 "I - 7* H - 1 

Proof of Lemma 7: On the one hand, if k n a r — > , , -, for all n G Z>n then 

xV = (k + a)"*' = E;=o ( ") fc"^' +r ^ E; = o ( •) (-l)"- j ^+T = on 

recognising the binomial expansion of (1 — 1)". On the other, if x n a r for all 

n e Z >0 then fc"a r = (x-a)"a r = E"=o (") * n ~ j ' (-l) J '-« j+r 4 (-l)"^+r 
on using the fact that only the j = n term is non-zero and invoking (145). Thus 
the result of lemma 7 holds in both directions. 

To complete the proof of lemma 6 we then deduce the pair of Cesaro limits 
for f n . r and g n ^ r given in lemma 7 simultaneously by induction on n. 

When ti = 0we have result (145) and when n = 1 we have 



P[f l,r](x) = 1 [ k+a /l,r(t)dt 
X Jo 




^ - (a; — a) (x — a — 1) + (x — a)a r+1 



x(r + 1) \2 
1 

2(7+1) 
1 



-{x-(2a + l) + 2a r+1 }+o(l) 



2(r + l) { r + 2) r + 2 

on invoking (145) and recalling that x is an eigenfunction of P with eigenvalue 
\ so that Climx = 0. This verifies lemma 6 and the formulae in (146) when 

x— >-oo 

n = 1 as our base case. 

Now, for the inductive step, suppose that gj r (x) — >• VI < j < n — 1 
and for all r > 0. Then in similar fashion to the working for the n = 1 case we 
have 
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k-1 



1 ™ * An 



. +1 r + 1 

3=0 



c(r + l) 

Now in the second term here, writing 

k n ■ a r+1 (x - a) n ■ a r+1 1 



--^ {b n+1 (k - 1) + k" ■ a^} (147) 



Z=0 

it follows by the inductive hypothesis that the only term with non-zero Cesaro 
limit is the I — 1 term, giving overall contribution to (147) of 

— ' ( U ) (-l) n_1 a n+r ~ (-l) n "S tt^ r (148) 

r + 1 V 1 / (r + l)(n + r + l) 1 ; 

For the first term, recall that b n +i{k — 1) = (B n +i (k) — B n+ \) and so 
we get overall contribution to (147) of 

I(r+1 ; (n+1) I' (T) 

If Z < n the inductive hypothesis applies to all except the j = 1 term, and when 
I = n + 1 it applies to all except the j = 1 and j = n + 1 terms. In the last 
case of i = j = n + 1, however, the resulting function is a pure power of x with 
generalised Cesaro limit in any case, and so the only contributions across the 
board come from the j = 1 terms, namely 



(r+ ^! l+1) E( n ! 1 ) B — (-iv 



(r + 1x1+1) {B n+ i(-l)-B n+ i} 
f-D" f-D™ 

:(« + !) = 7-"TTf ( 149 ) 



(r + l)(n+l) v y (r + 1) 
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Combining the contributions (148) & (149) in (147) it follows finally that 

P\f Ifr) 4 { - 1)n " n = { - ir 
[Jn ' riy ' (r + 1) (r + l)(n + r + l) (n + r + 1) 

and this completes the inductive step, thus completing the proof of lemma 6 
and the equivalent formulae in lemma 7. 
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